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ARE QUANTA UNIDIRECTIONAL? 
G. BREIT 


ABSTRACT 


Theory of the Brownian motion of a Planck radiator in black body 
radiation.—Jn the case of a classical resonator it is shown that the motion takes 
place even though the emission and absorption is not unidirectional provided 
the radiation is not strictly monochromatic. This theory is then extended to 
the case of a quantized atom. Interference between the waves emitted by the 
atom and the surrounding radiation provides for the recoil actions without 
the assumption of unidirectional quanta. Statistically these have the same 
effect as Einstein's unidirectional quanta. This theory removes the contra- 
diction postulated by Einstein between the Rubinowicz-Sommerfeld considera- 
tions and the requirements of classical Brownian motion. It leaves the 
wave-theory of light intact and is in agreement with Bohr’s principle of 
correspondence. 


INTRODUCTION 


XPERIMENTS on the development of photographic plates, on the 
photoelectric effect as well as the theoretical considerations of Ein- 
stein and Schroedinger! support the belief that quanta are darts of light, 
i.e. that the electromagnetic momentum of a quantum is in the same 
direction for all of its energy. Neither of these is a case of conclusive 
evidence. The photographic and photoelectric experiments indicate 
only the places where energy is absorbed and it appears that the only 
considerations dealing with the momentum of the quantum are the 
theoretical ones just cited. The fundamental point is contained in the 
paper of Einstein. Thisis verbally quoted asa proof of the unidirectional 
nature of quanta. In the paper there appears to be no statement to this 
effect though it shows that if it be granted that quanta are darts of energy 
1A, Einstein, Phys. Zeits. 18, 121, 1917 
E. Schroedinger, Phys. Zeits. 23, 301, 1922 
L. Silberstein and A. H. Trivelli, Phil. Mag. 44, 956, 1922 
L. Silberstein, Phil. Mag. July 1922, p. 257 


Svedberg and Andersson, Phot. Journ. August 1921, p. 325 
Svedberg, Phot. Journ. 1922, pp. 183-186 
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the conditions of correct Brownian movement are satisfied. The legiti- 
macy of the converse proposition is not evident—nor does the writer 
think it is true—though at first sight it is difficult to see in what way the 
conditions will be satisfied unless the hypothesis of light darts is made. 
In this note an alternative to Einstein’s treatment is presented which 
makes it unnecessary to assume the unidirectional nature of quanta. 


SIMILAR PROBLEM IN CLASSICAL THEORY 


This problem was first discussed by Einstein and Hopf? who showed 
that the requirement of correct Brownian movement of a classical 
resonator necessarily leads to the Rayleigh-Jeans formula. The discus- 
sion given below is not essentially different from that of Einstein and 
Hopf. It calls attention however to some features of the problem which 
are of importance in connection with the quantum case. 

Let us first see how classical electrodynamics explains the Brownian 
movement of a Planck resonator in a black body radiation. 

For simplicity let the resonator have an axis the orientation of which is 
fixed, and let its motion of translation be confined to that in a straight 
line which we choose as the axis of X (OX) and which we take as the 
resonator’s axis. We shall suppose the resonator attached to a particle 
of large mass M so that the translational velocity of the resonator is 
always small compared with that of light. The calculation will be similar 
to that used by Einstein. 

We consider the beginning and the end of a fixed time interval + and 
we change the beginning and the end at random keeping 7 constant. 
We take the mean square of the momentum for the beginning and for the 
end and we express the fact that the mean square is the same for both. 
If the velocity at the beginning of the interval is v, if the systematic re- 
sisting force is Rv and if the summation of the accidental random impulses 
is A 





(Mv — Rvr+A)? = (Mv)? 

and A?/r =2RkT (1) 

The origin of the resisting force Rv lies in the fact that there is a relative 
motion between the resonator and the radiation (the velocity v). The 
force Rv is of the same nature as radiation pressure and will have es- 
sentially the same expression in Einstein’s theory and this. 

The origin of A is not quite so simple. It lies in the spectral width of the 
band absorbed by the resonator because the absorption of one frequency 


? Einstein and Hopf, Ann. der Phys. 33, 1105, 1910. I am indebted to Professor 
Sommerfeld for this reference. 
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by a resonator put in the well known cube of Jeans’ does not give any 
recoil action on account of the stationary character of the waves. 

The resonator responds however not to a single frequency but to a 
band. Consider stationary waves of frequencies v;, ve The electro- 
magnetic momentum of the sum depends on [£,H,]+[E£.H.]+[£,4.]+ 
[Z.H,] where E, H are the electric and magnetic intensities corresponding 
to vy; and vg. The first two terms of the above expression correspond to 
the waves acting separately and give no flux during a time large in com- 
parison with the period. The last two terms represent the inte-ference 
effect of the waves. They give terms of the form sin (we—w,)¢ where 
w=2zv and thus correspond to long period oscillations in the electro- 
mag 1etic momentum. 

Calculation of A. For this we avail ourselves of a result of Planck.* 
Planck is concerned with an oscillator having a charge e, a mass m, exposed 
to radiation of energy density 67,dy in the frequency range (v,v-+dv). The 
energy of the oscillator itself is denoted by « and the amount of energy 
absorbed by the oscillator from the field is m,. Planck considers n? 


for a time 7 and takes its average value. The result is 
ne = —Tyer (2) 
2m 


Now the energy absorbed by an oscillator is the work done on it by the 


external field. In the time between two instants ¢;, 2 this work is 
ts 


= f cE,,fdt (3) 
h 
where E, is the component of the external electric field taken along the 
axis of the resonator and £ is the velocity of the oscillator’s electrified 
particle having the charge e. 

The force on the oscillator due to the action of the magnetic field on 
its moving electrified particle is perpendicular to the resonator’s axis 
and gives rise to the radiation pressure. The value of this force resolved 
along the line of motion of the resonator is (eH,/c)~ [OY being an axis 
perpendicular to OX and OZ] and using (2) and (3) 


Ba _Trer (4) 
2mc? 





since statistically the behaviour of E and H is the same (for justification 
see Note 1). 
Calculation of R. For this we use formula (260) of Planck*® 
ne=}(e?/m)T,r (5) 


3M. Planck, Wirmestrahlung, 4th edition p. 153, formula (261) 
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This gives the average absorbed energy in the time r. It may be shown 
that the radiation pressure is such as corresponds to the absorption of an 
electromagnetic momentum which is v’/c? times the absorbed energy 
where v’ (see Note 2) is v[3 —(v/p)(Ap/dv)]. In the case of the Rayleigh- 
Jeans radiation formula v’ =v and by (5) 


Ro = (eT,/4mc?)v (6) 
Using (4) and (6) it is seen that (1) is satisfied if 
e=kT (7) 


which is true in the classical theory. 

If the axis of the resonator forms an angle ¢ with the line of motion 
the terms A2/r, Rare found to be sin? ¢ times the values found above. 
Equation (1) is therefore again satisfied in virtue of (7). 

Resumé of treatment of classical case. The above derivation shows 
how a classical resonator reacting with the radiation by the ordinary laws 
of electromagnetism is maintained in a state of Brownian movement with 
the correct mean kinetic energy per degree of freedom (k7/2). The 
derivation is based on the similarity between the expressions for the 
energy absorbed and the radiation pressure. It also implies a statistical 
equivalence of EZ, and H, in the calculation of A. 

The above picture does not involve the hypothesis of unidirectional 
emission though it is true that this hypothesis also satisfies (1). It is 
clear that unless some additional fact is given the mere validity of (1) 
will not enable one to decide between these two possibilities and we are 
thus inclined to doubt the possibility of such a decision in quantum 
theory. 


DISCUSSION OF QUANTUM CASE 


Very little is known about the mechanism of absorption and emission 
of radiation. The discussion of the physical reality is of necessity only 
tentative. 

It will be satisfying to make the basic phenomenon of the quantum 
and the classical treatment the same. Interference between the field 
of the resonator and the surrounding field may be thought of as funda- 
mental in the classical theory and will be presented below as the basis 
of a quantum theory of Brownian movement. By so doing the apparent 
contradiction between the view of Rubinowicz and Sommerfeld, that the 
wave emitted by a Bohr atom may be spherical, and Einstein’s postulate 
of correct Brownian movement, is removed because the spherical wave 
hypothesis may be conceived as satisfying Einstein’s postulate in a man- 
ner analogous to that of the classical resonator. 
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For simplicity let us deal with an atom emitting linearly polarized 
light, (e.g. a quantized Planck resonator). It is feasible to suppose that 
the recoil action on the atom due to its emission or absorption is the same 
as if the same energy changes took place from a classical resonator. 

We can now calculate the terms A and R. For the latter the calcula- 
tion of Einstein can be interpreted to mean that v’ =v[3—(v/p)(dp/dv)] 
(see Note 2). Applying this to Einstein’s model and preserving the rela- 
tion c2A? =, 2 (see (2), (4)) it is found that both A?/7 and R can be obtained 
from his values by multiplying them by the factor 3. The relation (1) 
is therefore still satisfied. 

The fundamental postulates of this theory are then: 

(1) The systematic resisting force Rv has the same ratio to the absorbed 
energy as if the quantized atom were a classical resonator moving in the actual 
radiation field, here determined by Planck’s formula. 

(2) The relation c2A?=n2 is true not only in the classical but also in the 
quantum theory. 

In the above an atom emitting linearly polarized light, say a quantized 
Planck resonator was discussed. By interchanging electric and magnetic 
singularities the result is extended to circularly polarized light. 


The writer has spoken about this subject to several people whose ideas 
influenced his own. He first erroneously believed that one could think 
consistently of a quantized atom not exchanging momentum with the 
radiation. In a conversation Prof. Einstein was kind enough to point 
out that a classical resonator experiences radiation pressure and that it 
would be therefore inconsistent not to introduce an equivalent of this 
in the quantum theory. Thus it appears that the classical analogy has 
been kept in mind by Einstein but that he preferred to build the theory 
along revolutionary lines rather than to modify the classical theory to 
the least extent. The subject was also discussed with Prof. Ehrenfest of 
Leyden, Prof. Ames of Johns Hopkins, Prof. Epstein of Pasadena, Dr. 
Kemble and Dr. Van Vleck of Harvard. 


Note 1. STATISTICAL EQUIVALENCE OF E, AND Hy IN THE 
CALCULATION OF A?/r 


Since E, affects ¢ and H, does not, a justification is needed for the inter- 
change of these in { E,édt in the derivation of A2/r. 

The classical theory of black body radiation may be developed by 
considering the resonator as damped by its radiation resistance. To 
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within a good approximation this damping may be taken into account 
by a frictional term in its equation of motion‘ which we write 
E+KE+wo't = (e/m)E, 


xt is here the frictional term. Let 


°E,= f A (w) cos (wt — gy)dw. 
m 9 








é -(4 (w) cos {(wt— go—tan"{(w*—ao")/Kw)]} 5 
' V 2+ [(w? — coo?) /wo]? 


We are interested in the average square of the integral from #; to te of 
(e/m)E,fdt. Replacing integrations as to w by summations, multiplying 
the expressions and integrating with respect to the time, we obtain terms 
of the types A(w:)? and A(w;)A(w2). The latter type is subject to 
double summation as to w; and w2 and it contains a set of terms having 
in the denominator w;—w:2. When the result is squared these terms give 
the main contribution to the average if r=¢.—¢,; is a short time in com- 
parison with that necessary for the development of E, into the Fourier 
series. These terms give 





ye (w1)? A(we)? cos? ¥ sin? (w1—we)rt+4 sin? yp sin*t $(w; —we)r 
4 (wi — we)? {2+[(or2— 2/a1) FP} 

where y= — 9(w) — tan™{(w;? — wo?) /w1] + ¢(we) and the summation is taken 
over all values of w;, we excluding w;=ws. It is clear now that if A (we) be re- 
placed by B(we) and if A(we)?=B(we)? the result of the summation is 
unchanged on account of the random nature of ¢(ws) g(w:). This and a 
random change in the ¢’s is all that is required in order to replace E, 
by H,. Thus the validity of the interchange is justified. Taking the 
average and replacing the sum by an integral the average value becomes 
(x?/8)(A?/dw)?(7/x) where dw is the difference between two successive 
w’s in the summation. This is readily seen to be identical with the expres- 
sion of Planck. . 

The occurrence of the terms w;—w2 shows the importance of the inter- 
ference between adjacent frequencies for (A?/r). 





Note 2. THE CALCULATION OF THE RESISTING FORCE ON A CLASSICAL 
RESONATOR MOVING IN BLACK Bopy RADIATION 


Let E,’, Hy’ be the components of the electric intensity along the 
resonator’s axis and of the magnetic intensity perpendicular to the 


* Kretschmann, Ann. der Phys. 65, p. 310, 1921 
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resonator’s axis and to the axis of motion—both measured in a frame of 
reference moving with the resonator. The time integral of the radiation 
pressure is the time integral from #, to ¢# of (e/c)Hy,'é, where é is the 
velocity of the oscillatory particle. The absorbed energy is the corre- 
sponding integral of e&E,’'é. If the average over a long interval of time is 
desired the various frequency terms can be considered independently. 
The ratio of two terms of corresponding frequencies in both expressions 


ts 


is ai (1/c) Hy'(w)-E,'(w)dt)/( fE,'(w)dt) where the sign (w) after a 


h 

symbol indicates that only the term in the particular frequency 
discussed is taken and where it is supposed that H,’(w), E,’(w) 
arein phase. [The field is analyzed into a number of plane waves of 
different frequencies.| This ratio is the ratio of three times the 
density of electromagnetic momentum, resolved in the negative direction 
of OX, to that of the electromagnetic energy. Since resonance is sharp 
only a narrow band of frequencies need be considered. For this the 
ratio is the same throughout. 

In calculating the ratio we may consider a number of independent plane 
waves within the band sweeping by the resonator. No interference action 
between these waves must be considered because their frequencies are 
affected by Doppler’s principle to a different extent on account of the 
difference in direction. Thus we may use Eq. (19) of Einstein’ which 
gives for the energy in the frequency range dv per unit solid angle in a 
direction making an angle ¢ with the line of motion 


3 
~(6 + =. cos o7*) (: ~ cos ¢) ds 
4a c ~~ Ov c 


where pdy is the energy density with respect to the container. Multi- 
plying this by —cos ¢g/c, integrating, and dividing by }$pdv’ the ratio 
of radiation pressure to energy absorbed per unit time is 

(v/c*) [3—(»/p)(8p/ pr)] =v'/c* 
where v’=v[3—(v/p)(8p/dv)]. In the classical theory p is of the form 
av’? and v’ =». 


NATIONAL RESEARCH FELLOWSHIP, 
Crurt LABORATORY, HARVARD UNIVERSITY, 
March 1, 1923. 
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THE EFFECT OF A MAGNETIC FIELD ON THE 
ABSORPTION OF X-RAYS 


By JosepH A. BECKER' 


ABSTRACT 


Effect of a magnetic field on the absorption of x-rays.—A continuation of 
previous work (Phys. Rev. 20, p. 134, 1922) with improved specimen holders 
and technique. Jron, aluminum, and carbon showed no changes in the absorp- 
tion coefficient as great as one part in 10,000 (due to a magnetic field of 17,000 
gauss) when the tube was operated at 90 kv peak. Carbon, wood, and lithium 
also showed effects smaller than the limit of accuracy of the apparatus—2 parts 
in 10,000 for 61 kv, 4 parts in 10,000 for 46 kv, and 10 parts in 10,000 for 30 kv. 
If the effect exists, then, it is very small. 


HIS article is a continuation of a previous article of the same title 

in the Physical Review for August, 1922, vol. 20, p. 134. Unless 
otherwise stated the procedure was the same as described previously. 

The specific objects of the present investigation were (1) to design and 

test a new specimen holder to be used when the magnetic field H was 


perpendicular to the x-ray beam X; (2) to determine how closely the 
method of determining the relative change in the absorption coefficient 
due to a magnetic field could be relied upon; and (3) to extend the work 
to longer waves for a few low-atomic weight elements. 

The new specimen holder. In the previous article doubt was expressed 
as to the efficacy of holder h in Fig. 1 of the previous article, in preventing 
slight displacements of the specimen. The new holder was built along 
the lines used in the construction of holder h:, shown in the same figure. 
This new holder, which was used in all cases when H was perpendicular to 
X, rigidly held the specimen immediately in front of a lead slit or opening 
and prevented any mechanical effect which might be caused by the forces 
on the pole-pieces or by the direct action of H on the specimen. 

Duplicability of results. To determine whether or not the small 
effects obtained were due to chance and how closely the results could 
be duplicated, a specimen of iron, one of aluminum, and one of carbon 
were tested. The following week another series of tests was performed 
with the same samples. The same procedure was again repeated the 
third week. In all these tests the peak voltage across the tube was 90 kv 
and the magnetic field of 17000 gauss was perpendicular to the x-rays. 


1 National Research Fellow in Physics 
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Fig. 1 is a typical graph of the observations. The rate of drift of 
the electroscope leaf in divisions per minute 6 is plotted against the 
time in minutes. It will be noticed that the magnetic field was applied 
for one-half minute during each minute. The double arrow on the 
right shows Aé, i.e. the change in 6, that was produced by changing the 
area of one beam by four parts in ten thousand. This Aé could be 
determined with a precision of about 10 per cent. 


TABLE | 
Results for a field of 17,000 gauss, perpendicular to the x-rays 











Trial Ab 1/Io Al/I fe | tole 
no. due to H (x10-*) | (10-4) (10-4) 
Iron 1 + 2 45 oa 4 — .4+1 
2 +1.4 - 4 5 — .5+1 
3 —1.5 ° —.5 + + .6+1 
Aluminum 1 3 55 5 4 — .85+1 
2 1.4 - ian 4 — .5+1 
3 0 - 0 4 — O+1 
Carbon 1 4.2 .60 .6 4 —1.2+1 
2 —1.2 - —.2 3 + .4+1 
3 —0.1 —.0 4 0+1 


























Table I sums up the results of this series of tests. AJ/I gives the 
proportionate change in the intensity of the beam which has passed 
through the magnetized specimen. This is obtained by dividing Aé 
due to H by the Aé for a change in area of 1 part in 1000. The maximum 
oscillation is obtained by taking the difference in 6 for the highest and 
lowest points for H on, and H off, respectively, and dividing the larger 
of these two oscillations by Aé for A/1000. An effect equal to about 
one-fourth of this maximum oscillation could have been detected. Ay/y, 
the proportionate change in the absorption coefficient, is obtained from 
I/I, and AI/I by means of a formula developed in the previous article. 
The figures in the last three columns are reported as so many parts in 
10,000. 

It will be noticed that AJ/J is never more than 15 per cent of the maxi- 
mum oscillation. It can conservatively be concluded that a magnetic field 
of 17,000 gauss perpendicular to an x-ray beam produced by 90 kv peak 
changes the absorption coefficients of iron, aluminum, and carbon, by less than 
1 part in 10,000. 


LONGER WAVE-LENGTH TESTS 


Carbon, wood, and lithium were examined with x-ray beams produced 
by 61, 46, and 30 kv peak with the results given in Table II. Figs. 2, 3, 
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and 4 show representative graphs for wood at the three voltages. Thirty 
kv was as low a voltage as could advantageously be used with the present 
apparatus. 


TABLE II 








T/T; Al/I , Au/u 
(X10-*) (<10-) 


8 —1.6+3 


— 416 
+ .7243 


—2.4+4 
— 328 
—6.4+2 
—6.5+4 
+ .84+7 
+1.9+3 
—2.4+3 
—3.1+4 


. . 46 
Material H in gauss dus to 
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In general the effects are smaller than could be detected with certainty. 
In all but the two cases marked with a star in Table II, AJ/J is less than 
one-fourth of the maximum oscillation. The two exceptions are plotted 
in Figs. 2 and 3. During the course of the present trials it was attempted 





Aluminum Gkv HiAX Wood C6LA-. AWX 


Fig.t. 


Wood 45hkv. HIX Wood B30kv. ANX 





S=Dritt in Divisions per Minute. 














< 70 Minutes — 
Figs. 1 to4. e=H off;x=H on 


to keep conditions as steady as was possible with the facilities described. 
The smallest effect which could have been detected with certainty and 
measured with a precision of approximately 25 per cent was 2 parts in 
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10,000 at 61 kv, 4 parts in 10,000 at 45 kv, and 10 parts in 10,000 at 30 kv. 
To measure smaller effects than these would require a complete remodel- 
ling of the apparatus, paying particular attention to (1) steady voltage 
across and steady current through the tube, (2) constancy of the source 
of potential for the plates of the electroscope, (3) symmetry of the paths 
traversed by the two beams, (4) stray effects of the x-rays and magnetic 
field, (5) mechanical jars, (6) uneven temperature effects, and (7) the 
specimen holders. 


SUMMARY 


The results of the two papers may be summed up in the statement that 
magnetic fields up to 18,000 gauss do not change the absorption coef- 
ficient for x-rays generated in a tungsten target Coolidge tube by 90, 60, 
45, and 30 kv by more than about 1 part in 10,000 for the higher voltages 
to 1 part in 1000 for the lower voltages, in some ferromagnetic, dia- 
magnetic, low atomic weight and high atomic weight elements. 


ROCKEFELLER HALL, 
CORNELL UNIVERSITY,” 
March 21, 1923. 


? Present address: Norman Bridge Laboratory, California Institute of Technology 
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THE OPACITY OF AN IONIZED GAS 
By JOHN Q. STEWART 


ABSTRACT 


Theory of the absorption of radiant energy by free electrons in an ionized 
gas.—The classical theory of the optical properties of metals, supposed to 
contain free electrons, is applied to an ionized gas. As a result of collisions 
with molecules, vibrational energy supplied to the electrons by the electro- 
magnetic field is transformed into energy of thermal agitation. The coefficient 
of absorption comes out proportional to \*p*, where \ and p are the wave- 
length and pressure. There is also a scattering effect proportional to p and 
independent of \, but this is relatively unimportant except at pressures of 
10-‘atm. or less. The equations developed are tentatively applied to a dis- 
cussion of (1) opacity of the vapors of exploded wires. However Anderson 
reports an_opacity 200 times greater than that roughly computed from theory, 
and varying inversely as \ or \*. Evidently the subject requires further study. 
(2) Opacity of the solar atmosphere. Employing Saha’s theory to calculate the 
ionization, coefficients are computed which indicate that electronic absorption 
may be an important cause of the opacity of the solar photosphere. It is 
concluded that light from regions where the pressure is greater than .01 atm, is 
cut off completely, so that all we see comes from a spherical shell of rarefied gas. 
(3) Opacity of giant stars. Computations made, assuming the relation between 
temperature and pressure given by Eddington’s theory of stellar constitution, 
are found in rough agreement with Eddington’s values of opacity. 





EARS ago J. J. Thomson showed that, owing to scattering of radia- 

tion, an atmosphere containing free electrons should possess marked 
general opacity; and he suggested that this cause of opacity might be 
effective in the solar atmosphere.! Now that Saha’s theory is available 
for calculation of the ionization, this suggestion can be quantitatively 
applied. 

By the process of scattering, radiation is diffused in direction but not 
lessened in amount nor much changed in frequency. The term absorp- 
tion should be reserved for those processes which result in the disappear- 
ance of radiant energy. The present paper points out that in an ionized 
gas collisions between free electrons and molecules theoretically result 
in such an absorption. These collisions transform into thermal energy 
of translation the vibrational energy supplied to the electron by the 
electromagnetic radiation field. Formulas describing this process have 
been worked out by Drude, Lorentz, and others in connection with the 





1 J. J. Thomson, Phil. Mag. 4, 253, 1902 
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free electron theory of the optical properties of metals.? In the present 
paper these formulas are applied to conditions in an ionized gas.* 


THEORY 


The standpoint is that of classical electromagnetic theory. How far 
this theory represents the facts is doubtful. For this and other reasons 
the results of this paper are to be regarded as only tentative. 

When plane waves of polarized light are advancing through an ionized 
gas which is not too dense, the following well-known equation applies 
statistically to the motion of each free electron in a layer parallel to the 
wave-front and thin with reference to the wave-length, 

at meee + 2mru=eE. (1) 

dt c dt? 

The mass of an electron is m, or 8.93 X 10°"; its charge is e, or 4.774 K 10°"; 
its radius is a, =2e?/3mc?, or 1.88X10-"; c is the velocity of light, or 
3.00X10'°. (In this paper c.g.s. and electrostatic units are employed 
throughout, except that pressures are expressed in atmospheres.) The 
magnitude of the alternating electric vector of the radiation field is 
E;t represents time; and u is that part of the component of the electron’s 
velocity in the direction of E which is due to the action of the radiation 
field. The number of collisions with atoms or ions per unit time per 
electron, due to the thermal agitation, is r. 

The electric field of the radiation exerts a force eE on each electron; 
the direction of this force lies in the plane of the wave-front. There is 
an additional force due to the polarization of the medium; but this is 
negligibly small unless the density of the ionized gas is very great or the 
frequency of the radiation low compared with the frequency of light. 
The magnetic field exerts a very small force (the radiation pressure) 
in the direction of propagation of the radiation; and this force also is 
negligible in the present discussion. 

Eq. (1) states that the impressed force eE is balanced by three re- 
actions. The existence of the first is obvious. The second is the well- 
known reaction on an electron of its own radiation, and may also be 
written —(2e?/3c*) (d’u/dé#). The third is statistical; its existence is 
indicated by the following argument. 

If a steady electrostatic field E’ acts on an ionized gas, each free 
electron is subject to an acceleration eZ’/m in the direction of the field. 
The effect of collisions between electrons will be neglected. Since the 

* For example, O. W. Richardson, Electron Theory (1916), pp. 410, 432 


5A preliminary account of the work has appeared: J. Q. Stewart, Nature 111, 186, 
1923. 
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mass of an electron is relatively small, the magnitude of its velocity in a 
given direction after a collision with an atom may be assumed unin- 
fluenced by its velocity before the collision. Between successive colli- 
sions, then, each electron moves on the average in the direction of the 
field the distance eE’/2mr’; so that the average velocity of drift due to 
the electrostatic field is u=eE’/2mr. Thus the average momentum 
of each electron in the direction of the field is mu=eE’/2r. Statistically 
the effect of collisions, therefore, can be interpreted as a frictional force on 
each electron of amount 2mru. Since the velocity of drift is small com- 
pared with the average velocity of thermal agitation of the electrons, r 
is independent of wu. 

When E£’ is equal to unity, the magnitude of the electric current in 
a column of unit cross-section, or the electrical conductivity, obviously is, 
if there are m free electrons per unit volume, 

o =ne"/2mr. (2) 
The current carried by the comparatively massive ions is negligible. 

Returning to (1), write E=E,cos2rvt, where v is the frequency 
of the radiation of wave-length \. The average rate s at which energy 
is diverted by each electron from the primary radiation is (1 /t)SneEudt. 
The intensity J of the radiation is given by cE,?/8x. The quantity a/d 
is so small that its square may be neglected in comparison with unity. 
If r/v likewise is small, as is true in the cases to be considered, it follows 
from (1) that 

s=(3/mr) (2x*a/A+r/v)anrl. (3) 

If r=0, as for an isolated electron, all of the diverted energy is re- 
radiated—“‘scattered’”—with frequency unchanged (to the first order 
of small quantities); and no energy of the radiation is transformed to 
heat. The energy scattered per unit time per electron is‘ 

$;=62a?I = (8re*/3mc')I = (6.7 X 10°") J, (4) 
when s, and J are expressed in c.g.s. units. 

When r is not zero, the average rate at which heat is developed per unit 
volume is, from (3), 3nrac?E,?/8x*v?. This may be written 40,Eo, where 
o,, thus defined, is the electrical conductivity of the ionized gas for alter- 
nating current of frequency v. Accordingly, provided r/» is small, and 
making use of (2), 

oy = ne! /4x?mr*o, (5) 
nearly. 

The physical significance of (3) is better brought out by writing it in 
the form, 

s/r=(42°av/dr)6+2¢. (6) 
‘L. Page, Astrophys. J. 52, 67, 1920 
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Obviously s/r is the energy diverted per electron per collision interval. 
The quantity ¢ is the maximum kinetic energy of vibration of an electron 
due to the radiation; its value is easily shown from (1), when r/yv is small, 
to be (3ad?/2rc)J, since the maximum velocity attained in the simple 
harmonic vibration is, nearly, eEo/2mmyv. According to (4) the amount 
of energy scattered per electron per vibration is 
(6xa?/v)I = (42°a/d) o. 
Thus the first term on the right hand side of (6) represents the energy 
scattered by an electron during the time between two successive colli- 
sions. (Collisions are assumed instantaneous). Evidently, then, (6) 
implies that at the average collision, twice the maximum vibrational 
kinetic energy of an electron in the radiation field is transformed into 
disorganized heat-molecular energy. The appropriate assumptions about 
the nature of collisions should immediately give this result, without going 
through (1) with its statistical term 2mru. Another supporting line of 
argument has been given by Lorentz, in connection with the effect of 
molecular collisions on bound electrons. He has shown’ theoretically 
that if an electron is given r random blows per unit time the resultant 
term due to such “impact-damping” in the equation corresponding to 
(1) is 2mru. 
Eq. (3) can also be written 
$= (8mre*/3mc*)I + (2e*rd?/ame*) I, (7) 
or, in c.g.s. units, 
s=(6.7 X10°*)7+(16.0 x 10°) rA7J. 
The first term, or s;, expresses the amount of energy scattered, the second 
term Sz, the amount of energy absorbed, per unit time per electron. 
Write P 
K=K,+K2=ns;/I+ns,/I. 
Then K, is that part of the opacity coefficient K which is due to scatter- 
ing, and K¢ is the part due to absorption. 

It may be shown that the intensity J of radiation in the direct trans- 
mitted beam emergent from a layer of ionized gas of thickness z, when the 
intensity in the incident beam is Jo, is given by 

I=Ige*, (8) 
where ¢ is the base of natural logarithms. By (3) the ratio between the 
energy scattered in unit time per electron and the energy absorbed is 
2x*av/dr, a ratio which, owing to the decrease in 7, increases with de- 
creasing pressure of the gas. If scattering is appreciable the total amount 


5H. A. Lorentz, Proc. Amsterdam Acad. 8, 591, 1905, and other papers in the same 
Proceedings. 
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of light transmitted is greater than that which gets through in the direct 
beam. 

Before (7) can be applied it is necessary to obtain an expression for r 
in terms of known quantities. As the kinetic theory of an ionized gas is 
an unexplored field, the following argument may be greatly in error. 
Assuming equipartition of energy between molecules and free electrons, 
the root-mean-square velocity of a free electron in a gas is given by 
/(3RT/m), where R is the gas constant per molecule, or 1.37 10°", 
and T is the absolute temperature. This reduces to (6.8 10°)./T. 
The length of the electronic free path is 1/7NA?, where N is the number 
of atoms in unit volume, whether ionized or not, and A is the radius of 
an atom (assumed independent of ionization, an assumption doubtless 
incorrect). In most cases of astrophysical interest the ionized gas is 
presumably monatomic. Thus (nearly enough), 

r=nNA?/(3RT/m). (9) 

Write n= Ni, so that 7 is a numerical factor expressing the average 
amount of ionization. Also write N=273 Nop/T(1+7), where p is the 
total gas pressure in atmospheres, including the pressure 7p/(1+7) 
of the free electrons, and No is the number (2.705 X 10°) of molecules per 
unit volume of gas when p=1 and 7 =273. 
Then, in the units adopted, since 27347Nov/(3R/m) =1.6 X 10”, 





r=(1.6X108)A?p/(1+7)/T; (10) 
o =e /227A*+/(3mRT) =601/A?V/T; (11) 
Kya 2V 3273 NieVR A*\*ip*, (12) 
cm3/2 T*2(1 +1)? 
= (6.9 X 10°) A*n*ip?/T8/2(1+41)?; 
. K,=82net/3m?ct = (4.9 K 10" )ip/T(1+2), (13) 
since 82(273)Noe*/3mct =4.9 x 1073. 


When the mass coefficient of opacity k is employed, K in the above 
equations is replaced by kp where p is the density. Thus the mass 
coefficient of absorption corresponding to Ky is, for a monatomic gas, 

ke =2e?\*ri/amc?§MypW =4.4\*ri/W (14) 
= (5.7 X 10*8)A*Nip/W(1+i)/T, 
where W is the atomic weight, and Mp is the mass of the unit of atomic 
weight, or 1.65 X 10°. 

In order to apply to physical problems certain of the above equations 
it is necessary to know the approximate value of A, the collision radius 
of atoms or ions of the gas. In view of the tentative nature of the theory 
presented in this paper the determination of A need not be precise. 
The collision radii of the molecules of the permanent gases at ordinary 
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temperatures are known to be of the order 10°*cm. The atomic volumes 
indicate that this is the order of magnitude of the radii of the uncharged 
atoms of metals. A crude theoretical formula for the collision radii, with 
regard to electrons, of metallic ions may be derived by the following 
method. 

The ions, each carrying a positive change of average magnitude ie, 
may be regarded as centers of electrostatic forces which conform to 
Couloumb’s inverse square law. Accordingly, the energy equation of 
an electron moving in a hyperbolic orbit with velocity u at distance x 
from an ion is, neglecting the forces due to other ions, 

smu? = mu" +ie?/x, 
where xu’ is the velocity at infinity. Assuming equipartition of energy 
between electrons and other molecules, mu’ =3RT, on the average. 
It can easily be shown that the distance of closest approach between 
electron and ion, when the former is deflected through an angle 20 by 
the encounter, is (2ie?/3RT) cot?@. Thus, if @ is taken as 7/4 for the 
average collision, 


















A =(1.1X10-%)2/T. (15) 
As the ion is a structure as well as a center of force this value of A is to 
be regarded only as a rough approximation. 







ASTROPHYSICAL APPLICATIONS 














The opacity of the vapor of an exploded wire. The experiments of 
Anderson® have shown that a layer a few centimeters thick of the vapor 
of an exploded iron wire is sufficiently opaque to result in the production 
of a continuous spectrum, and to prevent the transmission of light. 
In a letter to the writer Dr. Anderson states that his most recent experi- 
mental results indicate that the value of K for iron vapor at 20,000° 
absolute and a total pressure of about 2 atmospheres, is of the order 
unity (presumably for \ near 4X10-); and, further, that K varies in- 
versely as \ or i’. 

The observed spectrum indicates that the iron atoms are for the most 
part doubly ionized, so that the value of 7 is about 2. Under these 
circumstances, (13) gives K,;=3.3 10-7; (10), (11), (12), respectively, 
give r=(7.5X10*)A*, o=0.85/A*, and Ke=(3.5X10")A?. The small 
size of K, shows that scattering is not important here. 

When i=2 and 7T=20,000, (15) gives A=10-’ cm. Accordingly, 
¢=8 X10" e.s.u./cm*, a conductivity seven thousand times less than 
that of metallic copper; and K,=0.004, which is 250 times smaller than 




















6 J. A. Anderson, Astrophys. J. 51, 37, 1920; Proc. Nat. Acad. Sci. 8, 231, 1922 
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Anderson’s experimental value. If Anderson’s published estimate of 
the pressure, namely ~/(1+7)=20 atmospheres, is employed, the cal- 
culated value of Kz comes out about 4. Thus, much depends upon 
the accuracy of the determination of the pressure. Perhaps a more 
serious disagreement between the theoretical K, and the experimentally 
observed K is that the former varies as \? while the latter, it is thought, 
varies inversely as \. Obviously the subject wants further study. 
The value of the “radius” of the iron ion, especially, is also very un- 
certain. 

Good electrical conductors are good reflectors, as well as good absorbers 
and emitters of light. If the boundary of the vapor of an exploded wire 
were sufficiently sharply defined a considerable portion of incident light 
might be reflected. Thus the ‘Heaviside layer” of ionized gas in the 
upper atmosphere of the earth has been supposed by radio engineers to 
reflect Hertzian waves. 

Pressures in the solar atmosphere. The following cursory examination 
of conditions in the outer regions of the sun indicates that the thermal 
absorption described by (12) may be an important cause of the opacity 
of the sun’s photosphere. 

The temperature of the sun’s photosphere and lower reversing layer 
is not greatly different from 6000°. Substituting in (13) gives 

10°°/K,=12(1+2)/ip (16) 
for the distance in kilometers which would be required to reduce the 
intensity of the direct beam in the ratio of 1/e if scattering by free 
electrons were the only cause of opacity. 

Ionization in the photosphere may not be sufficiently intense to make 
(15) applicable; but for purposes of argument A may be assumed as 
5X10-8. Then, for \=6X10-, (12) gives 

10-°/K,=(75 X 10°) (1+-2)?/ip? (17) 
as the distance in kilometers which would be required to reduce the 
intensity of the direct beam in the ratio of 1/e if the thermal absorption 
by free electrons were the only cause of opacity. 

Saha’s well-known theory permits calculation of 7, and consequently 
of 10-°/K, and 10°/K, in terms of p and T. Table I presents such 
calculated results for the case of sodium vapor at 6000°. The second 
ionization of sodium is neglected and Russell’s values of 7 are used.’ 
Values of r, the number of collisions per second per electron, are included 
from (10). This table may be regarded as illustrating, on a simplified 
scale, conditions in the solar atmosphere. The actual conditions are, 





7H. N. Russell, Astrophy. J. 55, 134, 1922 
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of course, much more complex, and the free electrons which are present 
have been liberated from all sorts of ions. 


TABLE I 
The opacity of ionized sodium at 6000°K 











Pressure in Fraction Collisions Kilometers distance to 
atmospheres, of of per second reduce intensity to 1/eth 
sodium vapor ionization per electron Scattering Absorption 
and free 
electrons 
p z r 10-5/K, 10-5/K; 

















1 0.21 4.2X10" 69 0.0052 
10-1 0.56 3.310" 330 0.33 
10-? 0.90 2.7 10° 2,500 30. 

10-* 0.99 2.5 X 10° 24,000 3000. 

10-4 1.00 2.5X 10" 240,000 300000. 











Table I brings out the predominating influence at moderate and high 
pressures of thermal absorption by free electrons in producing general 
opacity, as compared with scattering by free electrons. The latter process 
is relatively important only at extremely low pressures. It is worth 
recalling that electromagnetic theory shows that general, as distinct 
from sharply selective, scattering by bound electrons is not enormously 
greater (and may be much less) than scattering by free electrons in an 
appreciably ionized gas. 

Since K, varies as A’p’, a reduction in the estimate of the radius of 
the sodium ion would increase in the same proportion the values of 
pressure in the table. Ina crystal of sodium chloride the shortest distance 
between sodium and chlorine ions is known to be 2.8X10-§ cm. It 
seems unlikely that the collision radius of the sodium ion should be much 
less than half this. If thermal absorption by free electrons in the solar 
atmosphere really is as intense as the numbers in the last column of Table 
I indicate, it is at once apparent that at a depth where the pressure is as 
great as 0.01 atmosphere the opacity is sufficient to cut off light from 
lower levels. This figure, then, is indicated as an upper limit to the 
pressure in the visible regions in the solar atmosphere. Photospheric 
absorption may be supposed to set in at, or above, this depth. 

When observations are made of the flash spectrum, the grazing line 
of sight runs for tens of thousands of kilometers through the solar 
atmosphere. The faintness of the continuous spectrum under these 
conditions indicates, in connection with Table I, a pressure of less than 
0.001 atmosphere at these levels. A detailed examination of the condi- 
tions is outside the province of the present-paper. Other lines of ‘evi- 
dence, for example the sharpness of the Fraunhofer lines, indicate that 
the pressures are very low. 
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The opacity deep in stars. In his well-known theoretical study of 
conditions in the interior of stars Eddington® has found that the mass 
absorption coefficient has a value of the order 20. Recently he has 
dealt with this matter by a kind of quantum theory.’ 

It is of interest to consider the application of (14) to conditions in 
the interior of stars. If the wave-length considered is that of maximum 
energy radiation, the value of \ to be substituted in (14) is \=0.29/T. 
Making use also of (15), 


ko = (5.8 X 10”)i8b/WT®*(1 +2). (18) 
By Eddington’s theory of stellar constitution, 
(1.0 x 10°)p = BaT*/3(1+8), (19) 


where a is Stefan’s constant, or 7.6X10-", and the factor 1.0X10® is 
approximately the number of c.g.s. units of pressure in one atmosphere. 
For purposes of illustration the value of 8 may be taken as 0.83, which 
makes (19) apply to a giant star of mass 1.5 times that of the sun. 
Substituting in (18), 

ke =70028/W(1+2)V/T. 

Since by definition there are 7 free electrons for every atom, W may be 
written as W’/(1+7), where W’ is the average molecular weight of the 
positive ions. Accordingly, 

ke =70078/W'/T = 3078//T, (20) 
if W’ is 23. 

At the center of the star, T is of the order 410°, and 7 may be taken 
as 12, say. Thus k at the center is given as 25. Near the surface, say, 
where T is 10,000 and 7 is 2, ke is given as 2.4 Thus the opacity described 
by (20) is of the order of magnitude of the Eddington opacity, and also 
is roughly consistent with Eddington’s original assumption that k does 
not vary with depth in the star. 

In conclusion, the tentative character of the results presented scarcely 
needs to be emphasized. The subject of the opacity of gases to radiation 
is a broad and important one, which can be approached from many 
angles. In dealing with it the possibilities of classical electromagnetic 
theory appear by no means to have been exhausted. 

I wish to thank Professor Russell for his very helpful interest in the 
preparation of this paper. 


PRINCETON UNIVERSITY OBSERVATORY, 
March 7, 1923.!° 


* A. S. Eddington, Astrophys: J. 48, 212, 1918 
*°A. S. Eddington, Monthly Notices, 83, pp. 32, 98, 1923 
1 Received April 16, 1923 
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ON THE MOTIONS OF ELECTRONS IN GASES 
By K. T. Compton 


ABSTRACT 


Motion of electrons through gases when electric field and pressure are 
such that collisions are elastic.—For electron speeds below a critical value 
characteristic of each gas, the collisions with molecules are probably as between 
elastic spheres, since neither ionization nor radiation is produced. The 
average fraction of its energy lost in elastic collision with a molecule of mass 
M, by an electron of mass m, is f=2(m/M)(1—@/U), where Q and U are the 
average kinetic energies of molecules and electrons. But while, if U>Q, the 
electrons lose energy at collisions, they gain energy from the field between 
collisions, so the kinetic energy tends to increase or decrease to a terminal 
energy. Assuming Langevin’s equation for mobility, we find U; (in volts) 
=.66 El/\/f where El is the product of field strength and mean free path. 
In the case of plane parallel electrodes the number of collisions per cm and 
also the ratio of the speed at any point to the terminal speed, are both found 
to be independent of the field strength. In the case of a straight filament 
cathode surrounded by a coaxial cylindrical anode, the electrons may 
acquire their maximum speed near the cathode and subsequently lose speed 
while approaching the anode, their speed being greater than the terminal speed 
for this region. Curves are given for various cases. Substituting the value of 
terminal speed in Langevin’s equation leads to the following expression for 
electron mobility: pw =0.81514/ (ge/m) /[$2+ (22/44 W2)A]h where the energy 
due to the field is eW=leE/1.506\/f. This expression is somewhat similar 
to that given by Loeb (corrected for numerical error) but leads to values 
0 to 10 per cent higher. When the collisions are inelastic the mobility is 
shown to be greater. 


Bee a certain degree of approximation, which. varies with conditions, 
collisions of electrons with gas molecules may be regarded as if be- 
tween elastic spheres, provided the relative velocities do not exceed 
a value characteristic of the particular gas. The motions of such elec- 
trons in gases may then be treated by the ordinary methods of the 
Kinetic Theory of Gases. For many practical purposes such a treat- 
ment is sufficiently accurate. In this paper are presented some char- 
acteristics of electron motions in a gas in an electric field, with some 
applications to experimental results and a discussion of the conditions 
in which the method is inapplicable. We shall consider principally 
cases in which the gas pressure is sufficiently large, in comparison with 
the electric field, to make the average speed gained between successive 
collisions small compared with the total speed. 

Consider first electrons each of charge e, mass m, speed v, moving 
with a mean free path / in a gas under the influence of an electric field E. 
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If the speed gained between successive collisions is small compared with 
the total speed 3, all directions of motion are statistically equally proba- 
ble after an impact and the average distance advanced during each 
interval between collisions is 

s=}al?=}(eE/m) (t/i)?, 
where a=eE/m is the acceleration and //iv is the mean free time t. 

The mean rate of advance of the electrons, which may be written as 
the mobility u multiplied by the electric field EZ, is 

wE =s/t=}eEl/mi. 
Hence the electron mobility is 
u=sel/mi. 

This simple relation is inaccurate for several reasons. The effect 
of a distribution of free paths and of speeds cannot be averaged simply by 
use of average values of paths and speeds, and the dynamical effect of 
the collisions is such as to make the average initial speed after a collision 
different from zero. Townsend! gives an approximate treatment of the 
kinematical correction, which is more thoroughly treated by Lorentz; 
Maxwell’ considered the dynamical correction but not the complete 
kinematical one, while Langevin‘ and Lenard® considered both. Lange- 
vin’s equation for u may be writteA 


pod VG) ab m+M 
8/2 mc M 


where c is the square root of the mean square speed and M is the mass of a 
molecule. A critical discussion of these equations is given by Mayer.® 
Boltzmann’ and Langevin showed that the electrons acquire a Max- 
well’s distribution of speeds about their mean velocity of drift, so that 
c=1.086 v. Since m is very small compared with M, we may write (1) 
u=0.75 el/mi, (2) 
in which the factor 0.75 differs from that quoted by Townsend, 0.815, 
because of our substitution of i for c. 
Pidduck® has recently shown that the electrons do not acquire exactly 
a Maxwell distribution of speeds, and that the departure from such a 
distribution may increase the numerical factor by about eleven per cent. 





(1) 


' Townsend, Electricity in Gases, p. 83 

* Lorentz, Theory of Electrons, p. 68 

* Maxwell, Scientific Papers, 1, 398, 1890 
* Langevin, Ann. de Chim. et de Phys., (8), 5, 245, 1905 
5Lenard, Ann. der Phys. 60, 329,1920 

® Mayer, Jahr. Radioakt. Elek. 18, 201, 1922 

7 Boltzmann, Gastheorie, vol. 1, p. 114 

8 Pidduck, Proc. London Math. Soc. 15, 89, 1916 
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In the following treatment, we shall neglect this departure and assume 
the correctness of equations (1) and (2). 

Since wE is the average distance advanced in the direction of the 
electric field in a second, and i/} is the average number of collisions made 
during the same time, 

s=lwE/v=0.75 PeE/m?* (3) 
is the average advance between successive collisions. If we put i= 
0.849 c?, as in a Maxwell’s distribution, and eU =4mc?, we obtain 

s=0.441 PE/U. (4) 
Thus the average number of collisions 1/s made while advancing unit 
distance varies inversely as /? and hence directly as the square of the gas 
pressure p. 

If vy is the number of electrons per unit volume, the current density 
j of electrons is 


j=vepE =0.75le’Ev/mi (5) 
and the density of negative space charge p is 
p=ve=jmv/0.75 leE. (6) 


NUMBER OF COLLISIONS MADE WHILE ADVANCING A 
DISTANCE FROM THE CATHODE 


Consider first electrons emitted with negligible speed from a plane 
cathode and drawn through a gas toward a plane parallel anode. The 
total average number of collisions per electron in distance d is 


N= -f' ten f Bete, Exdx 113. (7) 
0.441 PE p 


This differs by the hater 1.134 from a similar expression derived by 
Hertz.° 

If the electrons are emitted from a cylinder of radius r,, such as a wire 
filament, and are drawn toward a coaxial cylindrical anode, the average 
number of collisions made by an electron while attaining a distance 


from the axis is 
r r 


v= f U dr= 
4 0.441PE 0.4417 


since, if o is the charge per unit naa of the siiiiiite E=20/r and 
=f Edr=2clog r/r.. This integral gives 


N= 1.1347} og —a4y't. (8) 
P r r 


—dr 











* Hertz, Physics, 2, 15, 1922 
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This reduces, of course, to Eq. (7) if r=r.+d, where d is small compared 
with ¢,, i.e., when the distance from the cathode is so small compared 
with its radius that it can be considered as a plane surface. 

Eqs. (7) and (8) are inaccurate in that they assume no loss of energy 
by the electrons as they advance through the field. Such losses actually 
exist, for even if the collisions are perfectly elastic, some momentum is 
transferred to the molecules at each impact. At high gas pressures these 
energy losses, though individually very small, are in the aggregate quite 
large and limit the speed which the electrons attain in any given field. 
The consequences of these energy losses are considered in the following 
two sections. 


ENERGY LossEs AT COLLISIONS 


If an electron collides with a molecule, initially at rest, and is thereby 
deflected through an angle ¢ with a change in speed from v to 2, the 
molecule receives a velocity w; in 
the direction @. Assuming perfect 
elasticity (as for electron impacts 
below the critical speed for in- 
elastic impacts) we have the re- 
Fig. 1 lations 





mv-+mv, cos ¢= Mw, cos 6 
mv, sin ¢= Ma, sin 6 
mv? —4mv,? =} Mw;’, 
from which 

bmv?—pmv,? _ 4 Mm 
mv? (M+m)? 
gives the fraction of its energy lost by the electron at the collision. Actu- 
ally electrons may strike the molecule at any angle 6 between 0 and 7/2, 


and the average fraction of the energy lost is 
/2 


(a | t02s sin 00080 a= 2—____—__ 


cos’6 


fo= 





Mm 
(M+m)? 
Since M is much larger than m, we may put 

f=2m/M. (9) 

We have assumed, thus far, that the molecules are initially at rest. 
We may take into account the effect of initial molecular velocities by 
averaging the effects of two types of collisions, between electrons and 
molecules moving in opposite directions with average speeds c and w 
and between those moving in the same direction. Out of N collisions, 
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there are 3}N(c+w)/c of the first type, resulting in an average energy 
loss equal to M[(m/M)*c?+(m/M)cw], and }N(c—w)/c collisions of the 
second type, with an average energy loss equal to M[(m/M)*c?—(m/M)cw}. 
We can therefore obtain the total loss of energy by the electrons in 
all N collisions, and thus the average loss per collision. When this is 
divided by the average energy before collision we obtain.'° 

f=2 mS )=27 ae | (10) 

Me M 

If the speed c of the electrons is less than that which would make them 
in thermal equilibrium with the gas, f is negative and the electrons gain 
energy at encounters, on the average. If, as is the case in an electric 
field, their mean kinetic energy exceeds that of the gas molecules, they 
lose energy at collisions. Eq. (9) may be used with sufficient accuracy 
except under conditions, such as high gas pressures and small fields, in 
which the mean energy of the electrons is not much greater than that of 
the molecules. 

Eq. (9) has been experimentally verified for helium and it is shown in a 
later paper that it is probably true for impacts in nitrogen and hydrogen. 
We may assume it to be correct for impacts of electrons with monatomic 
molecules at speeds below their minimum radiating speeds, and for at 
least some multiatomic molecules under similar conditions. 


TERMINAL SPEEDS 


In advancing a distance dx in a field E each electron gains energy eEdx, 
and loses an average of feU at each of the dx/s intervening collisions. 
The net gain of energy per electron is 

ed U =e(E—fU/s)dx. 
Substituting for f and s from (10) and (4), we have 
dU/dx = E—4.536 mU(U—Q)/PME. (11) 
The terminal speed of steady drift is found by solving this equation 
for U when dU/dx =0, giving 
tntoe / FE ” 
4 4.536m 
in equivalent potential drop: In terms of average speed, since }mc?=eU 
and c= 1.086%, and the mean energy e2 of the molecules is given in terms 
of absolute temperature by a7, 


| ‘ee 
i= [ert (or + i) ] (13) 
1.086 m 1.134 m 








1° Benade and Compton, Phys. Rev. 9;+87,4948- 
4 34. 1913 
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In most practical cases, 2 or aT are so small compared with the other 
terms that they can be neglected. This is equivalent to using Eq. (9) 
in place of Eq. (10) for f. In such cases Eqs. (12) and (13) for terminal 


speed reduce to 
hese / M (14) 
2.13 m 


Wi 
v= 0.801 (%* 4/H) (15) 
m m 


We may obtain an approximate idea of the magnitudes of these 
terminal speeds under various conditions by putting the mean free path 
1=4,/2-L, where L is the mean free path of a gas molecule. This rela- 
tion is known to be approximately correct, although there are deviations 
from it, which will be discussed in a later paper. From Table I the ter- 
minal energies in equivalent volts may easily be computed for various 
fields and gas pressures in the case of some common gases. 


TABLE I 


Terminal energies in various gases 


1, =electronic free path at 1 mm pressure and 20°C. 
U:=A E/p, where E is in volts/cm and p in mm. 








Gas /(M/m) 1,(cm) U, (volts) 
He 85.9 1313 5.30 E/p 
A 271.5 -0461 5.86 
H2 60.8 .0842, 2.40 
N2 227.4 .043 4.65 
CO; 285.0 .0290 3.88 
Hg 608.0 0135 3.85 








These, values | of U; are meaningless (1) if E and p are so chosen as to 
make Ue comparable with the mean energy of thermal agitation Q= 
0.0372 “volts, i in “which case Eq. (12) must be used in place of Eq. (14); 
(2) if E’ and” _p are so chosen as to make U; comparable with the critical 
energy “at which electron impacts become inelastic. Between these 
values Table I should be approximately correct. 

When in this terminal state of speed, the average number of collisions 
made per centimeter advance is given by Eqs. (4) and (14) as 

_ 1.005, /M 





(16) 


MEAN ENERGY OF ELECTRONS AT ANY DISTANCE FROM CATHODE 
(a) Plane parallel electrodes. The rate of gain of energy at any point 
is given by Eq. (11). Putting 
4.536 m/PM =4.536 p*m/1l,°M =a? (17) 




















MOTIONS OF ELECTRONS IN GASES 339 


and integrating from U=0 at the cathode x=0 to U=U at the distance 
x=x from the cathode, we find the mean energy at distance x to be te ahah 


u=yat ey 148 tev (1/e+ Ms /4E ~1< Tha ‘ea 
es 6 41 
If the sz E= 0, this reduces to U=Q, i.e., the energy of the electrons 
is simply that of thermal agitation. If the energy due to the field is so 
large that 2 may be neglected, Eq. (18) becomes 
_Eeé*—-1 
a +41 
which is identical, except for a small correction in the value of a, with 
an equation previously derived." When x is very large, U approaches 
the terminal value U;=E/a, which is identical with Eq. (14). 

The mean energies of the electrons at any distance from the cathode 
for various values of field and gas pressure with different gases are shown 
in Figs. 2 and 3._ Fig. 2 gives U as a function of x for various fields and 
values of a, and Fig. 3 indicates the type of gas and its pressure corre- 
sponding to any value of a. 

The average number of collisions made while going a distance x from 
the cathode is given by (4) as 


w= ftar- f° as 
y 8 : 0.441 PE 


If we substitute for U from Eq. (18) and integrate we obtain the correct 
number of collisions. If we substitute from Eq. (19) we get a number 
which is quite accurate except in very weak fields. From Eq. (19) we 
find readily that 
N=(M/4m) [log (2+4%*+ 6?) —log 4]. (20) 
This equation reduces to N=0 if x=0. If x is very large N approaches 
x times the value of 1/s given by Eq. (16) for the case of terminal speed. 
If energy losses are negligible, which means M/m has a very large value, 
N takes the value given in Eq. (7). 
From Eq. (19) we can readily compute the distance d in which the 
electrons acquire a fraction ¢ of their terminal energy to be 
se ite (21) 
2a 1-¢ 
and the average number of collisions made while acquiring this fraction 
of the terminal energy is 


(18) 





(19) 








N= oe : 
4m i-,¢’ 


(22) 


1! Benade and Compton, loc. cit.’ 
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Similarly, the average time to go a given distance or acquire a given 
fraction of the terminal energy may easily be computed from these equa- 
tions. It is interesting that the number of collisions made in any given 
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22 44 
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1.6 32 
mA 28 
2 
> 24 
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Fig. 2. Energy of electrons as a function of distance from a plane cathode for certain 
values of the field E and pressure factor a. 
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Fig. 3. Pressures corresponding to various values of a for various gases. 


distance is independent of the strength of field; the distance in which a 
fraction ¢g of the terminal energy is acquired is independent of the field; 
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and the number of collisions made while acquiring any fraction of the 
terminal energy is independent of both the field strength and the mean 
free path. 

The terminal speeds are acquired surprisingly quickly in gases at con- 
siderable pressure, as illustrated in Table IT. 


TABLE II 
Distance in which fraction ¢ of terminal energy is acquired (0°C) 











Helium Nitrogen Mercury 
p ¢ d (cm) d (cm) d (cm) 
760 0.1 0.00064 0.00057 0.00047 
0.2 0.00131 0.00115 0.00096 
0.5 0.00356 0.00313 0.00262 
0.9 0.00953 0.00840 0.00708 
0.99 0.01710 0.01507 0.01260 
10 0.1 0.049 0.043 0.036 
0.9 0.724 0.638 0.533 
1 0.1 0.49 0.43 0.360 
0.9 7.24 6.38 5.33 








(b) Coaxial cylindrical electrodes. This case is of particular interest 
because it corresponds to the common arrangement of a hot wire fila- 
ment source of electrons situated at the axis of a surrounding cylindrical 
anode. - 
The rate of gain of energy by an electron at any point between the 
electrodes is given by Eq. (11), with 
E=2¢/r and 20=(V,—V,)/log (ra/r:), (23) 
where V,, V, and r,, r, are the potentials and radii of the anode and 
cathode, respectively. Introducing a from Eq. (17), we have 
dU _ 2 a°QUr—a?U*r 
dr r a 2o 

I have been unable to find a solution 

of this equation, but have evaluated 

the integral by a laborious step-by- 


step process, as follows. Ku 





(24) 


If r; and re are very nearly equal, 
we may consider the field in the Soe 
interval rz—r;, to be constant and 
given by E=+/(E£,£:;) and calculated 
from Eqs. (23) in terms of the applied 
potential difference and the dimen- 
sions of the electrodes. We may 
then use 
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dU/dr =(E*+a°QU —a°U")/E 


and integrate to find the increase in U through this range from 1; to r2 
Us Lé 


f dU a f dr 
/ E+aQU-aU? J E 
I have solved the special case where the term in Q is negligible, as is the 
case when electron energies are larger than about 0.05 of the minimum 
ionizing energy of any gas. We then find the mean energy U; at the end 
of the interval in terms of the energy U, at the beginning of the interval 
to be 
















U,= E(E+aU,)e*"-") — (E—aU;) 
= 

a(E +a U)e20"—") + (E +a U,) 
This process may be carried out over each of the small intervals into which 


, (25) 





a . @& E 8 9 
T (cm) 


Fig. 5. Energy of electrons at various points between cylindrical electrodes. 





the space between the electrodes is subdivided, and the energy at every 
point thus determined. 

Examples of the way in which electrons gain energy between cylindrical 
electrodes in gases are given in Fig. 5. The gas pressure for any type of 
gas corresponding to each value of a may be found in Fig. 3. The 
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curves are drawn for cases of 4, 10 and 20 mil wires (r.=0.005, 0.0125, 
0.0250 cm, respectively) with an anode of 1 cm radius and an applied 
potential difference V,— V,=20.66 volts. The same curves may be 
made to apply to other cases by altering the units according to the rela- 
tions shown in Eq. (23). 

It is seen that even at quite moderate gas pressures, as low as 1 mm, 
the aggregate loss of energy at impacts is quite large and the electrons 
fail to acquire the energy corresponding to the applied potential differ- 
ence. At the higher pressures the electrons acquire their greatest energy 
quite close to the filament, and lose energy as they proceed toward the 
anode into the regions of less electric‘intensity. The maximum value of 
U for each curve is that which is equal to the terminal speed for the 
electric intensity at that particular point, as given by U=E/a, by Eq. 
(14). To the left of this maximum the electrons have less energy than 
the terminal value corresponding to their location, while to the right 
they have greater energy than the corresponding terminal value. This 
is illustrated by the curve of theoretical terminal energy drawn for one 
particular case (a=4). 

These facts have a bearing on the design of ionization tubes for work 
at low voltages. For instance, the theory of cumulative ionization™ 
indicates that the effect of imprisoned resonance radiation in promoting 
ionization should increase with increase of gas pressure and thus an 
increase of gas pressure should facilitate the striking of a low voltage 
arc. This is true up to about 5 mm (depending on the electrode distance), 
but above this pressure the arc strikes with increasing difficulty and 
higher voltage is required. This is to be expected in view of the large 
energy loss by impacts at the higher pressures. 

Since, at these gas pressures, the electrons have approximately a Max- 
well’s distribution of speeds, some may ionize even though the mean 
energy is below the minimum ionizing energy. The greater part of the 
ionization will occur in the region of the maximum value of U. 

The above considerations still hold in principle, but the equations and 
numerical values must be altered, if the electron current density is so 
great as to appreciably alter the potential distribution between the 
electrodes. In the case of a large current from a filament, for instance, 
the negative space charge shifts the regions of large electric intensity 
toward the anode, so that curves like those in Fig. 5 would reach a lower 
maximum value of U and this maximum would be shifted toward the 
anode. The correct curves giving the mean energy at every point 


® Compton, Phys. Rev. 20, 283, 1922; Phil. Mag. 45, 1923 
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could be obtained by the method of Eq. (25) if the actual distribution 
of potential between the electrodes were known. 


ELECTRON MOBILITIES IN GASES 


The electron mobility in electric fields is given by Eq. (2) provided 
the mean speed 3d is so small that there are no appreciable number of 
ionizing or otherwise inelastic impacts. (According to Pidduck® the 
numerical factor may be as large as 0.84.) Substitution for the average 
speed from Eq. (13) gives, for the mobility, 

_ 0.815el 
© VmlaT + (o?T?-+PME%/1.134m) 
If the speed due to the field is negligible compared with that due to ther- 
mal agitation, this equation reduces to the familiar mobility equation™ 
wo =0.815 el/+/(2maT) =0.815 eb /mco=0.75 el/mio. (27) 
If, on the other hand, the ratio E/p is so large that the speed of thermal 
agitation is negligible compared with that due to the field, 


y= 0.842 / —% _ 9,707 4/ 2¥F (28) 
EV/(Mm) Em 


When E/> is small, the mobility is therefore independent of the field; 
when E/p is large, the mobility decreases as the field increases. The 
mobility is evidently greater the less elastic are the collisions, as is shown 
by the presence of the quantity f in the numerator of Eq. (28). In terms 
of f, the complete mobility Eq. (26) takes the form 
7 0.815 el 

eV maT + (@2T?-+1.76 PE%/fyy 
The so-called ‘‘mobility constant” K, is defined as the mobility of 

electrons in the gas at 760 mm pressure and 0°C, and is given by 
K =(p/760) (273/T)4u. (30) 
Loeb" was the first to derive a mobility equation in which the effect 
of both thermal agitation and the electric field were considered in es- 
timating the mean velocity of drift to be substituted in Eq. (2). His 
equation, when expressed in the symbols of the present paper, is 

a 0.815 el 
*  Vm|2aT +0.575 lEe/ Vfl 

The difference between Loeb’s equation and Eq. (29) may be seen most 
easily by putting a7 =eQ and JEe/1.506./f=eW, where e@ is the energy 
of thermal agitation and eW is the energy which the electrons would have 





(26) 





(29) 





(31) 


18 Townsend, Phil. Mag. 40, 505, 1920 
4 Loeb, Phys. Rev. 19, 24, 1922 
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acquired from the field had the thermal energy of the molecules not been 
contributed. In this notation Eqs. (29) and (31) take the forms 
isile 0.815 el 
V (2me) (42+ (2/44) 
a 0.815 el 
r /(2me)[Q2+0.433 W}' 

The presence of the factor 0.433 in Loeb’s equation is due (1) to his 
use of a value of the terminal energy which I derived in an earlier paper” 
and which is less accurate than the value used in the present paper and 
(2) to his use of this energy as a maximum rather than a mean value, 
thereby erroneously dividing it by 2. 

Apart from this numerical inaccuracy, which is accidental, the two 
expressions differ in form. Loeb takes the resultant terminal energy of 
electrons to be the sum of the energies which would be due to thermal 
agitation and to the field if each agency contributed additively. Thus 
his term in brackets, after correcting the numerical error, is U’=2+W 
whereas the correct value is 

U=430+(9°/44+ 0")! (32) 
If we put k equal to the ratio of the actual energy U to the energy of 
thermal agitation 2, we have from (32) 





(29a) 





(31a) 


U=kQ, W/2Q=~+/[k(k—1)], (33) 
whereas on Loeb’s assumption we would have 
U'=k’Q, W/Q=k’—1, (34) 


which is the same as an equation derived as an approximation by Pid- 
duck and frequently used by Townsend.*® 

It is found that the equations of Loeb and Pidduck are correct at 
both extremes, when k= 1 and when k=, but that there is a small error 
at intermediate values, reaching its maximum at k=4/3, when W= 
(2/3). Here k’ is 25 per cent larger than the true value k, and the 
mobility u’ calculated by Loeb’s corrected equation is about 10 per cent 
less than the theoretical value yw. It is likely that some of the conclu- 
sions regarding small departures from perfect elasticity of electron 
impacts in gases, which Townsend has based on Eq. (34), should be 
reconsidered in the light of Eq. (33). 


SPECIAL CASE OF COMPLETELY INELASTIC IMPACTS 


This case cannot be considered as a special case of elastic impacts, 
using Eq. (31) with f=1, since that equation was derived on the assump- 


% Benade and Compton, loc. cit.!° 

% It should be noted that we have used the symbols W and © as energies, whereas 
Pidduck and Townsend use them as speeds, so that they appear squared in their equation 
corresponding to (34). 
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tion that the speed gained between impacts is small compared with the 
mean speed. In the case of perfectly inelastic impacts we assume that 
the electron starts from rest after each encounter. The mean speed of 
advance is easily found as follows. 

Consider the N electrons which cross the plane x=0 in any long 
interval of time. Ndx/I of these collided previously in the layer between 
x and x+dx. Of these, the fraction «*/’ has arrived at the plane with- 
out making an intervening collision. Thus (N/I)e*/dx electrons arrive 
directly from distances between x and x+dx. Their kinetic energy is 
3mv?=Eex, whence x=mv?/2Ee and dx=(mv/Ee)dv. Substituting 
above, we have 

F(v)dv = (mv/lEe)e "7 = dy (35) 
as the fraction which arrives with speeds between v and v+dv. The 
average speed is 


oo 


i= fs F(v) dv= 4/ Ee. (36) 
V 


2m™ 


»=4/f .. (37) 
2 mE 


It is readily seen that the mobility in these circumstances is consider- 
ably greater than it would be if the impacts were elastic. Probably 
impacts in all gases are elastic below critical velocities characteristic of 
each gas. These critical velocities are known to be much higher than 
the mean thermal velocity in most cases, although there are probably 
some gases in which electron impacts are inelastic at speeds below ordi- 
nary speeds of thermal agitation. 

A test and discussion of these mobility equations will be presented 
shortly in another paper. 


PALMER PuysiIcAL LABORATORY, 
PRINCETON, NEW JERSEY, 
May 1, 1923. 





The mobility is »= v/E, or 


















CURRENTS BETWEEN COAXIAL CYLINDERS 


CURRENTS LIMITED BY SPACE CHARGE BETWEEN 
COAXIAL CYLINDERS 


By IrvinG LANGMUIR AND KATHARINE B. BLODGETT 


ABSTRACT 


Limiting current between coaxial cylinders; calculation of the function 
8=f(r/r.) in the space charge equation i =(2./2/9) /(e/m) V?/2/(rB?). Two 
different infinite series were obtained for 8 and the coefficients for fourteen 
terms of each were determined. These two sets were checked against each 
other. Thus values of 6? were obtained accurate to 1 in 10,000. These were 
checked by an integration method which was also used to calculate values in 
the region where the series failed. For a cathode of radius ro inside an anode 
of radius 7, 8? increases from 1 at r/ro=11.2 to a maximum value 1.0946 at 
r/r9=42, decreases to a minimum value 0.9990 at r/r).=30,000, and becomes 
1 at r/ro=0. The customary assumption §*=1 has therefore led to errors 
up to 9.5 percent in previous calculations, but this error is usually about 
equal and opposite to that introduced by neglecting the effect of initial veloci- 
ties. For the cathode outside the anode, 8? is given very closely by the equa- 
tion 8?=4.6712 (ro/r) [logio(ro/r) —logiox/2]*/2, for ro/r>10. The empirical 
constant 4/2 is interpreted to mean that the potential distribution near the 
anode is unaltered if the hot cathode is replaced by a cold cylinder having 
one half the cross-section of the original cathode. The correction for initial 
velocities is less for a cylindrical cathode inside an anode than for parallel 
planes. In the inverted.case it is much greater than in the case for parallel 
planes, and the effect of the tangential component of the initial velocity may 
greatly decrease the current that flows. 





N PREVIOUS papers! it has been shown that in a very high vacuum 
the space charge of the electrons limits the current that can flow 
from a hot cathode to a positively charged anode so that beyond a certain 
point an increase in the temperature of the cathode causes no further 
increase in current. When the cathode and anode are long coaxial 
cylinders the current limited by space charge is given by 
j2V2,/@ V8 
9 m (rB?) 
Here 7 is the electron current per unit length along the axis, V the voltage 
at any point P, r the radius at P, e and m the charge and mass respectively 
of an electron. 
Expressing 7, 7, and V in amperes per cm, cm, and volts, respectively, 
the equation becomes 


(1) 





i= 14.68 X 10-°V2/ (rB?) (1a) 


1 Langmuir, Phys. Rev. 2, 450, 1913; Phys. Zeits. 15, 348, 1914 
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8=f(r/ro) is a quantity of zero dimensions, ro being the radius of the 
cathode. Values for 6? were given in a short table, but these were derived 
from a rough method of calculation and many of them have an error of 
several per cent. The object of the present paper is to correct and extend 
this table and to give values of 8 for the inverted case where the cathode 
is the outside cylinder. 

It was shown in the derivation of Eq. (1) that 


r(@V/dr*) +dV/dr=i - (2) 
e 
Substitution of Eq. (1) in Eq. (2) gives 
36r?(d°B/dr*) +r°(dB8/dr)?+76r(d8/dr) +6*—1=0 * (3) 
Placing 
x =log (r/ro) (4) 
we obtain 
38(d°B/d-y*) + (dB/dy)?+48(dB/dy) +6°—1=0 (S) 
The solution of this equation gives the series? 
B=y—(2/5)y*+(11/120) 7° — (47/3300)yY'+ .... (6) 


This was the series used to calculate 8? but is accurate only for small 
values of y. The values obtained are all too low, the error at r/ro=5 
amounting to —.011. 

The value 8 = 1 is a solution of Eq. (3) and corresponds to the case that 
r/ro=. A-second method? for calculating 8 was based on a substitution 
B=1-—e* 
and a series was obtained for ». The values of 8 obtained in this way 
converged so rapidly to unity for r/r9=11 to 20 that the accuracy of the 
method was held to be doubtful. This series was therefore abandoned in 
favor of a method of extending the series in Eq. (6) to give accurate values 
of 8 for large values of y. 


DERIVATION OF SERIES FOR 8 





The four terms of the series in Eq. (6) fail to give values of 8 accurate 
to five places of decimals at r/ro>=1.7.. A much improved method sug- 
gested by H. M. Mott-Smith has made it possible to determine ten more 
terms of this series. This method consists in expanding 8, which was 
given as a function of y in Eq. (5), in the form 


B=AgtAwytAoy’+Asyt+ .... tAnv*t.... (7) 








2 Due toa printer’s error the coefficient of the second term was incorrectly printed 
—%% in the Phys. Rev. paper, but it was correctly given in a paper in the Phys. Zeits. 
J. A. Fleming pointed out this error in the Radio Review 2, 133, 1921. 

3 Langmuir, Phys. Zeits. 15, 351, 1914 
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From Maclaurin’s series we have 


= Betin( ® (Se of “(=) a 
B= Bot( =) +5 =3),t +25) + (8) 


where the terms with the subscript 0 have values corresponding to y =0. 
The values of the successive derivatives can be conveniently calculated 
as follows: 
Let B=a, dB/dy=b, d*B/dy*=c, etc. 
Then Eq. (5) takes the form 
3ac+b?+4ab+a?—1=0 (9) 
- Differentiating (9) we have 
3ad +5bc+4b? +4ac + 2ab =0 (10) 
Differentiating (10) we have 
3ae+8bd +-4ad +-12bc +5c*?+2ac+2b?=0 (11) 
Differentiating (11) we have 
3af+11be+4ae+18cd + 16bd +2ad+12c*?+6bc =0 (12) 
From Eq. (6) we see that 8=0 when y=0. Therefore in the present 
expansion for 8 in the neighborhood of y =0 we may place a9=0 and Eq. 
(9) reduces to 
be—1=0, do=1 (13) 
the positive sign being chosen since 8 increases with y and hence the first 
derivative is positive. Substituting a9=0, bo =1 in Eq. (10) we have 
Co= —4/5 (14) 
Substituting the values of ao, bo, and co in Eq. (11) we have 
dy= 11/20 (15) 
Similarly from Eq. (12) we have 
éo9= —376/1100 (16) 
By comparing (7) and (8) it is seen that 


Ao =) (17) 
n!\dy"/o 


Therefore dividing the derived values of bo, co, dp and eo by 1, 2, 6 and 24 
respectively we obtain 

A,=1, Ap=—2/5, As=11/120, Ag= —47/3300 
and these coefficients agree with those in Eq. (6). 

It is seen that in this process when each successive equation is dif- 
ferentiated the new equation thus formed can be solved for a value of 
(d"B/dy")9 by a knowledge of the values of all the preceding derivatives. 
Thus the coefficients of any number of terms of the series in Eq. (7) can 
be calculated. The 2nd column of Table I gives the first fourteen of these 
coefficients calculated by this method. 
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H. M. Mott-Smith also suggested a second series, usually more rapidly 
convergent than the first, which he obtained by putting 


B=067, (18) 
which reduced Eq. (5) to 
36 “+(S) -e=0. (19) 
dy? \dy 
This equation was also expanded by Maclaurin’s series, using the same 
method as that described above. 


TABLE [| 


Coefficients of terms in series for B 








B= ZAny" B =e-Y/2>Bny” 
An Bn 


= 





.0 
0 
| 


{0:10 

.09166667 + .01666667 
.01424242 + .002424242 
.001679275 + .0002872294 
.0001612219 + .00002658476 
.00001293486 +1.766124 x 10-* 
. 87693 X 10-7 +6.332946 x 10-8 
.46192 x 10-8 —8.73852 X 10°" 
-94843 X 10-° — 1.93844 X 107"! 
. 36026 x 10°” +5.77287 X10“ 
.1101 X10-" +9.4502 x10-" 
.6644 10°" +4.7012 X10-" 
.2526 X10-§ —6.5539 X10" 


CaNOU WN © 








When combined with Eq. (18) it takes the form 
»B=eV?(BotBytBayt+ .... +Biy"+...) (20) 


1 n 
em (*) (21) 
n!\dy"/Jo 


The third column of Table I gives the first fourteen coefficients and 
shows that this series which is at first more convergent than the A 
series becomes less convergent after the eleventh term. 

Each coefficient of the series was checked against the corresponding 
coefficient of the A series by the following method. Differentiating Eq. 
(18) successively and placing y=0 and hence 8B=0 and 6=0 we obtain 
a series of equations for d"8/dy" in terms of d@/dy, d?0/dy?, . . . . d"6/dy". 
From Egs. (21) and (17), substituting~ !B,=d"@/dy" for each deriva- 
tive of 0, and n!A,=d"B/dy" we obtain A,=B,; 2A.=2B.—B,; 6A;= 
6B;—3B,+%B;; etc. Each coefficient was tested by means of these equa- 
tions using one or two more decimal places than are given in Table I, 


where 
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and all were shown to be correct within the degree of accuracy recorded 
in the table. 


CALCULATION OF 8 


Case where r/ro>1. (Cathode inside of anode). The B series was used 
for all calculations of the value of 8 when r/ro>1 because it is the more 
convergent series over the first part of the range, while over the latter 
part it has the advantage that the error of the neglected final terms is 
multiplied by the increasingly small quantity «7/? which becomes 0.1 at 
r/ro=100. The fourteen terms were sufficient for an accurate calcula- 
tion of the values of 8 up to y=4.2, r/ro=66.7, and were useful as an 
approximation as far as y=7.5, r/ro= 1808. 

In order to check the results obtained, the following method was 
adopted. For y=3.0 to y=4.2, 8 was calculated for equal intervals 
Ay =0.2, and in addition the value of d8/dy was calculated at each point 
from the equation obtained by differentiating (20). The values of 8 and 
dB/dy for each value of y were then substituted in (5) and each equation 
solved for d?8/dy?. The successive values of d?8/dy? thus obtained were 
integrated by Simpson’s rule to give d8/dy and again integrated to give B, 
the integration in each case being checked by Weddle’s rule. It was 
found that the values of 8 obtained at the end of this process as a result 
of integration were exactly the same as those substituted in the equation 
at the beginning of the process. It was therefore shown that the series 
calculations for 8 (requiring all fourteen terms at y=4.2) had yielded 
correct values which were solutions of the differential equation for 8. 

When the series began to fail at y =4.2 the integration process proved 
to be a practicable method of determining the correct values. For if 
preliminary values of 8 and dB/dy, which are approximately correct, be 
substituted in the equation to obtain d?8/dy? and the integration process 
carried through as before, the result of integration will be a new set of 
values of d8/dy and 8 slightly different from the first. These new 
values may in turn be substituted in the equation and the process 
repeated, and if the final values of 6 do not then agree with the values 
obtained in the preceding approximation the process may be carried 
through as many times as necessary. It is found that the values of 8 
thus successively obtained converge toward a fixed set of values which 
are the desired solutions of the equation, so that the method provides a 
means of arriving at these solutions. 

From 7 =4.5 to y=7.5 the series calculations were used as preliminary 
values for the integration method, and the results of the operation of the 
method showed that these values of 8 required at y =4.5 a correction of 
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—0.00003, and at y=7.5 a correction of —0.00260. Beyond 7=7.5 
the series was abandoned and the preliminary values obtained by extra- 
polation. 

The results of these calculations for y=3.0 to y=10.5 are given in 
Table II, together with the corresponding values of d8/dy. 


TABLE II 








d3/dy 


+0.03506 
.02239 
.01221 
.00415 
.00211 
.00686 
.01035 
.01369 
.01527 
.01562 
-01509 
-01402 
.01263 
.01107 
.00949 
.00797 
-00656 
-00530 
-00352 
.00221 
.00126 
.00061 
.00020 
.00001 
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Table III gives the complete range of values of 6?, which were all calcu- 
lated to at least one more decimal place than appears in the table. The 
values corresponding to r/ro>20 were derived from interpolations for 8 
from Table II, using Newton’s interpolation formula. They show that 
8? passes through unity at r/ro=11.2, increases to a maximum value 
1.0946 at r/ro=42, and again passes through unity close to r/r9= 10,000. 
A minimum is reached just before y=10.5 (r/ro=36,316) and its value 
differs from unity by less than 0.0012. Probably from this point onward 
the curve passes successively through an infinite number of maxima 
and minima which rapidly approach unity. 

It has been customary to place 6?=1 in the space charge equation 
(1a), but the data in Table III indicate that in the range of diameters of 
cathode and anode most commonly used in electron tubes (r/ro between 
20 and 500) this practice has involved an error up to about 9.5 per cent. 
The reasons that this error has not been detected by experiments are 
probably (1) that the ‘end convections” are usually not accurately 
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known and (2) that this error is usually about equal and opposite to that 
introduced by neglecting the effect of initial velocities.‘ 


TABLE III 


8° as function of radius 
ro =radius of cathode; r=radius at any point P 
B? applies to case where P is outside cathode, r> 1. 
(—8)? applies to case where P is inside cathode, ro>r. 








r/t9 or oof r B (-8)? r/o or fo/r 


0.00000 0.00000 
.00010 
.00039 .00040 
.00159 
.00356 
.00630 
.00980 
.02186 
.03849 
.08504 
. 14856 
‘ . 2282 

1525 . 3233 

. 1854 .4332 

.2177 
. 2491 
2793 
. 3083 
3361 
. 3626 
. 3879 
4121 
.4351 
.4571 
-4780 
.4980 
.5170 
.5526 
.5851 
.6148 
.6420 
.6671 
.6902 
7115 
7313 
. 7496 
. 7666 
. 7825 
.7973 
.8111 
.8241 
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1.2 
1.3 
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1.5 
1.6 
1.7 
1.8 
1.9 
2.0 
2.1 
2.2 
2.3 
2.4 
2.9 
2.6 
2.7 
2.8 
2.9 
3.0 
3.2 
3.4 
3.6 
3.8 
4.0 
4.2 
4.4 
4.6 
4.8 
5.0 
5.2 
5.4 
5.6 
5.8 








Case where r/ro<1. (Cathode outside of anode). For the inverted 
case of electron current between concentric cylinders where the cathode 
is the outside cylinder, the same series may be used as in the preceding 
case, but y is now negative so that alternate terms of the series have 


- *This has been discussed in a footnote to a recent paper, Langmuir, Phys. Rev. 
21, 435, (1923) and illustrated by Dushman’s experimental data. 
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changed sign. This makes the values of 6 negative, but since only 6? 
occurs in Eq. (1) the current remains positive. 

The calculations of 8 were made from the B series as far as ro/r = 20. 
From this point they were carried through to the end, using the A 
series, for after the factor €’/? in the B series reached the value /20 it 
multiplied the error of neglected terms by too large a quantity. Also 
the A series is more convergent in this region. The results are given in 
Table III and are believed to be accurate to the number of places given. 

For high values of ro/r a useful approximation formula for obtaining 
6? is given by 

6?=4.6712 (ro/r) [logio (ro/r) —0. 1505]*/” (22) 

This is accurate to one part in 10,000 for values of ro/r greater than 
80, and accurate within one per cent for values of ro/r greater than 10. 
The equation is obtained from Eq. (19) by neglecting ¢’, integrating 
twice, substituting the value of @ into Eq. (18) and eliminating y by Eq. 
(4). The two integration constants are then so determined as to make 
the values of 6? for large ro/r agree with those calculated from the series 
as given in Table III. 

Logarithmic differentiation of Eq. (1), considering 7 as constant, gives 

3(r/V)/(dV/dr) =2+4(d log B)/(d log r) (23) 
Logarithmic differentiation of Eq. (22) gives 
4(d log 8)/(d log r) = —2—3/[log,(ro/r) —0. 34654] 
Combining this with Eq. (23) we have 
dV/dr = —(V/r)/[log.(ro/r) —0. 34654] (24) 
This equation for the potential gradient at any point is of course valid 
only in the same range as Eq. (22), that is, when r is less than say 0.170. 

In the absence of electron current and the resulting space charge the 

field intensity between concentric cylinders is 

dV /dr =(V/r)/log(ra/r) (25) 
where 7, is the radius of one cylinder at zero potential and V is the poten- 
tial at any point of radius r. 

Comparing Eqs. (24) and (25) we see that the field intensities become 
identical for all values of r if we place 

log (ro/r,) =0. 34654 =logV/2 
or ro=teV2 (26) 

The ratio ro:7, is equal to \/2 within about one part in 30,000, and 
therefore probably is rigorously equal to tite number. This result may 
be interpreted as follows. 

With current limited by space charge between an outer hot cathode 
and an inner concentric anode of small diameter, the potential distribu- 
tion near the anode is unaltered if the cathode is replaced by a cold 
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cylinder having one half the cross section of the original cathode. The 
effect of space charge is thus to increase the field intensity near the 
anode to the same extent as a 50 per cent reduction in the cross section of 
the outer cylinder in the case where space charge is absent. 

From this it appears that the space charge close to the anode is without 
appreciable effect in limiting the current. 

Effect of initial velocities of the electrons. In the derivation of Eq. (2) 
it was assumed that the electrons leave the cathode without initial 
velocities. In a previous paper’ a full discussion has been given of the 
effects produced by these initial velocities in the case of parallel plane 
electrodes. 

With a cylindrical cathode of small diameter in a large coaxial anode 
there is a relatively high potential gradient near the cathode so that the 
electrons move with high velocity through the large part of their paths. 
The correction for the initial velocities is thus considerably less than for 
parallel planes (about one fourth). 

In the inverted case, with the cathode as outer cylinder, the electrons 
move with low velocity over most of their paths. The correction for 
initial radial velocities must therefore be much greater than for parallel 
planes. 

In the inverted case the tangential velocities of the emitted electrons 
must be of very great importance in determining the current limitation 
by space charge. Consider an electron which is emitted from the cathode 
with a tangential velocity corresponding to a fall through the potential 
Vo. Under the influence of the central force the electron describes an 
orbit in which the angular momentum is constant. The tangential 
velocity thus increases in proportion to the reciprocal of the radius. 
Since voltages are proportional to the squares of the corresponding 
velocities, we have. 

Vi=Vo(ro/r)? 
where V, is the voltage equivalent to the tangential velocity at a point 
of radius r. 

If V; is the voltage of the anode with respect to the cathode, it is clear 
that (V,+Vo) cannot exceed V;. Neglecting Vo in comparison with 
V,, we see that the only electrons which can reach the anode are those for 
which 

Vo< Vi(r/ro)? (27) 

Taking the energy corresponding to the tangential velocity as kT /2 

and placing this equal to Voe we find 
Vo=kT/(2e) =7/23,200 (volts). 


5 Langmuir, Phys. Rev. 21, 419, 1923 
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With adsorbed films of alkali metals® it is possible to obtain large elec- 
tron currents at temperatures as low as 700°K under conditions particu- 
larly favorable for the use of a cathode outside of an anode. It is thus 
possible to make Vo (for the average tangential velocity) as low as about 
.03 volt. With this value of V> and the anode voltage V; equal to 300, 
we see by Eq. (27) that ro/r<100. Thus if the diameter of the anode is 
less than 1/100th of that of the coaxial cathode, about half the emitted 
electrons cannot reach the anode because of their tangential velocities. 
Since these, however, contribute about equally to the space charge, the 
electron current flowing to the anode will be approximately half that 
calculated by Eq. (1). These considerations make it possible to estimate 
the magnitude of deviations from Eq. (1) resulting from tangential veloci- 
ties and will serve as the basis for future theoretical and experimental 
investigations. 

Lack of radial symmetry in the electrodes is another factor that can 
cause tangential velocity components and lead to a decrease in the space 
charge current. 


RESEARCH LABORATORY, 
GENERAL ELECcTRIC Co., 
SCHENECTADY, N. Y., 
March 14, 1923. 


* Langmuir and Kingdon, Science, 57, 58 (1923) 
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ABSTRACT 


Electron emission from tungsten filaments containing thoria.—Ten years 
ago the erratic behavior of some tungsten filaments was traced to the thoria 
present and it was discovered that by suitable treatment the filaments, con- 
taining 1 to 2 per cent of thoria, could be activated so as to give an electron 
emission many thousand times that of a pure tungsten filament at the same 
temperature. Extensive series of experiments have led to the conclusions 
summarized in this paper. The increased emission is due to a layer of Th atoms 
adsorbed on the surface of the tungsten. To form this layer, the thorium 
oxide which originally exists throughout the volume of the filament, must first 
be reduced in part, by heating for a few seconds to between 2600° and 2800° K; 
then, by heating for some time at a suitable activating temperature 7, 
=2000° to 2100° K, some of the Th atoms are brought to the surface by 
diffusing faster than they evaporate, the fraction of the surface covered 
6, increasing to a limit @.« which may or may not be 1. The activity deter- 
mined at a lower testing temperature 7;=1400° to 1500° K, increases loga- 
rithmically with 6; for any filament it remains constant provided the filament 
is not heated above 1900° K and the film is not allowed to be oxidized. It is 
evident that the rate of activation d@/dt depends directly on the cifference 
between the rate of diffusion to the surface and the rate of evaporation from 
the surface. Rate of diffusion of Th atoms to the surface is DG where G is the 
concentration gradient (atoms per cm‘) of Th atoms near the surface, and D 
is the coefficient of diffusion. If No is the number of Th atoms per unit area 
for a saturated film, it seems probable, from the crystal structures of Th and 
W, that No is only half the number of tungsten atoms in the surface; and 
making correction for the increased surface due to the dodecahedral structure 
(6 per cent), No is taken as 0.75610". A mathematical analysis of the 
experimental results gives for the absolute value of D (cm? per sec.): logisD = 
.044 —20540/T, in agreement with a theoretical expression given by Dushman 
and Langmuir. Rate of normal evaporation of Th atoms from a partly covered 
tungsten surface (@=.2 to .8) in atoms per cm? per sec. is given by: logioEam = 
31.43—44500/T. As the surface covered decreases below .2, E, decreases 
but not as fast as @, while for values of 6 above .8, E, increases toa value three 
times Enm. The small variation of E, with 6 shows that the effect of adjacent 
Th atoms is small. Studies of the effect of prolonged activation of filaments 
indicate that the surface film is never more than one’atom thick. Since the 
atoms do not diffuse back into the tungsten, we must suppose that when a Th 
atom diffusing to the surface arrives below another Th atom, one Th atom 
immediately evaporates. The rate of induced evaporation depends upon the 
diffusion and upon 6, experiments showing that E;=DG(0.820+.186). 
The rate of activation is, then, given by: Nodé/dt=DG—E,—E;= 
DG(1—0.820—.18@)—E,. The experimental curves, however, correspond al- 
most equally well to the simple expression: d6/di =k(@«—6), when @ is the 
limiting value, at which the diffusion and evaporation are balanced. Deactévation 
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occurs when the evaporation exceeds the diffusion. The concentration gradient 
G depends, of course, on the previous history of the filament, especially the 
preliminary reduction of thoria at a high temperature. The rate of reduction 
of thoria is given by: logiop; =27.98 —30160/T, where , is the number of 
atoms of thorium produced per sec. per cm’; the limiting concentration gradient 
(atoms per cm‘) reached is given by: logiwG,=25.22 —9620/T, the actual 
concentrations of thorium in parts per million by weight varying from 1.4 to 
200 for temperatures of 1800° to 3000° K. Various properties of thoriated 
filaments in a steady state are summarized in Table VIII. 

Energy changes for atoms of thorium, per gram atom, computed from 
constants in the above equations are: (1) Heat of reduction of thoria in tungsten 
= —138,000 cal. (absorbed); (2) Heat of evaporation of Th from tungsten 
= 204,000 cal.; (3) Heat of diffusion of Th through tungsten: 94,000 cal. 

Crystal structure of tungsten and chromium.—Both are body-centered 
cubic and when etched show rhombic dodecahedron faces (110). Both show 
pronounced cubic cleavage (100). 


ARLY in 1913 Mr. William Rogers and the writer were endeavoring 
to determine the electron emission from tungsten filaments in the 
highest possible vacuum. For this purpose two hairpin shaped filaments 
were mounted side by side in a tubular bulb. Two platinum wires were 
sealed through the glass walls on opposite sides of the bulb, the wires 
being bent so that the part in the bulb lay flush with the inner surface 
of the glass. After baking out the bulb for several hours at about 400°C 
and glowing the filaments, the bulb was sealed off, and one of the filaments 
was then heated for a considerable time at about 2900°K until enough 
tungsten had evaporated and deposited on the inner surface of the bulb 
to reduce the resistance measured between the platinum wires to about 
100 ohms. Thus a gas free anode was formed and residual gases were 
removed by the clean-up effect of the vaporized tungsten. 

The electron emissions obtained in this way were considerably larger 
than those previously found under less perfect vacuum conditions. 
Since earlier work had shown that the clean up is more effective at low 
bulb temperatures, the bulb was now immersed in liquid air and more 
tungsten vaporized. With this still better vacuum very great increases 
in electron emission were obtained, often of many thousand fold. The 
results even in this vacuum were often erratic and sometimes changed 
rapidly with time. 

It was then found that these increases in emission depended upon the 
presence of small amounts of thorium in the filament. In making fila- 
ments for tungsten lamps it was customary at that time to add a fraction 
of one per cent of thorium oxide to prevent the occurrence of a phenome- 
non known as offsetting. Measurements with specially prepared pure 
tungsten filaments showed that the emission remained unchanged when 
the bulb was immersed in liquid air, and was, in fact, the same as when 





EMISSION FROM THORIATED FILAMENTS 359 


the emission was measured in presence of argon at voltages high enough to 
ionize the gas. The lower values of the emission from the thoriated fila- 
ment, obtained for example before the improvement of the vacuum, 
agreed fairly well with those from the pure tungsten. 

During 1913 and 1914 and again especially during the last four years, 
very detailed studies of the mechanism of the emission from the thoriated 
filament have been carried out.' In the following pages the main conclu- 
sions reached in this work will be given, but it will only be possible to give 
a small fraction of the experimental evidence upon which those are based. 

Conditions for obtaining high emission from a thoriated filament. If a 
thoriated filament is used in place of a pure tungsten filament in ordinary 
electron tubes, the chances are that nothing unusual will be noted in 
measuring electron emission, for the thorium on the surface evaporates 
rapidly at the normal operating temperature of a tungsten filament and is 
very sensitive to contamination by gases which may oxidize it. The 
thorium may be protected by using an exceptionally high vacuum such 
as that obtained by immersing the bulb in liquid air, or by vaporizing into 
the bulb such substances as sodium, potassium, magnesium or calcium, 
or traces of vapors of reducing carbon compounds such as hydrocarbons. 
Magnesium or calcium are most easily introduced by placing small 
pieces of these metals on the anode or cathode where they may be 
vaporized during the operation of the tube, or they may be placed on (or 
in) a small helical tungsten filament mounted in the bulb. Experiments 
have shown that the emission of a fully activated thoriated filament is 
the same whether the filament is in a very high vacuum or in sodium or 
similar vapor. 

In the ordinary thoriated tungsten filament nearly all the thorium is 
present in the form of thorium oxide ThO,. Apparently for this reason 
high emissions are ordinarily not obtained from thoriated filaments until 
after these have been heated to a very high temperature such as 2600°K 
ormore. A flashing of the filament at 2800°K for 3 minutes is a sufficient 
preliminary treatment. This probably reduces some of the thorium oxide 
to metallic thorium which remains distributed at very low concentration 
throughout the filament. The heating to such a high temperature, 


1 The effect of thorium in increasing the electron emission from tungsten was de- 
scribed by the writer before the American Physical Society on October 31, 1914 (See 
Phys. Rev. 4, 544, 1914) but was not published. A full description of the conditions 
needed to obtain this effect, together with an explanation of the phenomena, was 
given in U. S. Patents 1244216 and 1244217 issued to the writer on October 23, 1917. 
Very brief descriptions of some of the observed effects were published in Trans. Amer. 
Electrochem. Soc. 29, 137, 1916; Jour. Amer. Chem. Soc. 38, 2280, 1916, and Phys. 
Rev. 20, 107, 1922. 
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however, removes all thorium from the surface of the filament by evapora- 

tion, so that if the temperature is lowered to say 1900°K or less and the 
emission is measured it is found to be the same as from a pure tungsten 
filament at the same temperature. 

To activate the filament metallic thorium must be brought to the surface 
by diffusion from the interior. This occurs at an appreciable rate only 
at temperatures above about 1900°K. On the other hand at temperatures 
above 2200°K thorium evaporates rather rapidly from the surface of the 
filament. The best temperatures for activating the filament are there- 
fore in the neighborhood of 2000°K or 2100°K. 

After the filament has been brought to a very active state by this 
thermal treatment it will remain indefinitely in this condition if the 
temperature is kept below about 2100°K, provided the vacuum is suf- 
ficiently high or oxidizing gases are eliminated by the presence of easily 
oxidizable substances such as alkali vapors. In presence of inert gases, 
etc., the cathode may be deactivated by positive ion bombardment if the 
positive ions are given high velocities produced for example by high 
anode voltages.?, With lower voltages, the thoriated filament may be 
kept in an active state in presence of large pressures of argon, alkali 
vapor, etc., even when the gas is ionized sufficiently to eliminate practi- 
cally all space charge. 

Measurement of the activity of the surface. The activity of the filament 
under various conditions is studied by measuring the electron emission 
(saturation current) at some fixed temperature below that at which 
appreciable diffusion or evaporation occurs. We shall call such a tem- 
perature the testing temperature and represent it by 7;. If it is desired to 
study states of activity ranging all the way from pure tungsten to a fully 
activated surface, a testing temperature of 1500°K is most convenient, 
for at this temperature the emission from pure tungsten is about 0.1 
microampere per cm? while that from the fully activated surface is 18 
milliamperes per cm’, and this is a good range of current for use with a 
galvanometer. 

The laws governing the activation of a thoriated filament are best 
determined by means of activation curves such as that illustrated by the 
curve OAB in Fig. 1. The typical experimental procedure for obtaining 
these data is as follows. 

Three filaments a, b, and c, each of about 0.0075 cm diameter and 9 cm 
long, were bent into horseshoe shape and mounted side by side upon 
a single glass stem by means of platinum leads. This stem was then 
sealed into a tubular lead glass bulb about 13 cm long and 4.5 cm in 


* See K. H. Kingdon and I. Langmuir, Phys. Rev. 20, 108 (1922) 
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diameter provided with two flush platinum seals at opposite sides of the 
bulb near the centers of the filament loops. The glass stem on which the 
filaments were mounted had an enlargement blown in it between the seal 
and the flange which by casting a shadow prevented tungsten from de- 
positing on the surface of the flange, thus avoiding leakage currents be- 
tween the cathode filaments and the tungsten deposit which served as 
anode. 

The bulb was heated several hours to 360°C while being exhausted by a 
glass condensation pump through a trap immersed in liquid air. During 
this time the filaments were heated for short times to 2000°K. After 
cooling the bulb to room temperature the filaments were aged for a few 
minutes at 2500°K. Each of the filaments was then made anode while 
charged to +260 volts, the other two filaments being heated to 2500° 
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Fig. 1. Rates of activation and deactivation of a thoriated tungsten filament. 


and serving as cathode. This treatment heated the platinum leads by 
electron bombardment. The tube was then sealed off, the bulb was 
immersed in cold water (or often liquid air) and the central filament (0) 
was heated at 2900° until tungsten had deposited on the walls in such 
amount that the resistance between the platinum flush seals was only 
a few hundred ohms. This tungsten deposit, connecting the two platinum 
seals in parallel, served as anode in all measurements of electron emission 
from any filament used as cathode, the unused filaments being connected 
to the negative end of the cathode. For these measurements the anode 
was kept at 250 volts, this being sufficient to give saturation with the 
currents used. The bulb was always kept completely immersed in liquid 
air when the electron emission was measured. 

The filament to be studied, which usually contained from 1 to 2 per 
cent of thorium oxide, was flashed at 2800° for 30 seconds or more, and 
was then ready for obtaining data for an activation curve. After such 
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measurements the filament needed to be deactivated before another 
activation curve could be obtained. For this purpose the filament was 
heated for a certain time tz to a deactivating temperature Tg at which thor- 
ium evaporated more rapidly from the surface than it diffused to it. 
Temperatures from 2200° to 2600°K were usually chosen as deactivating 
temperatures. 

To determine the rate of activation of the filament at a given tempera- 
ture the filament was first deactivated and then heated for successive 
periods of time at the activating temperature 7,. After each of these 
periods the activity of the filament was determined by measuring the 
electron emission at the testing temperature 7;. By plotting the elec- 
tron emission i; measured at the testing temperature 7; (say 1500°K) 
as ordinate, against the time ¢, during which the filament has been kept 
at the activating temperature (2100°K) as abscissa, a curve is obtained 
like that shown by OAB in Fig. 1. The current at first is so small that 
it cannot be shown in a curve such as Fig. 1 although by using the more 
sensitive galvanometer shunts the currents are readily measured. After 
a certain time the curve rises very rapidly about like an exponential 
function, but then rises more slowly and finally approaches asymptoti- 
cally a limiting value. 

At any given point, such as A in Fig. 1, the activation of the filament 
may be interrupted and a deactivation curve AC may be obtained by 
heating the filament for successive short intervals of time at a deactiv- 
ating temperature 7, which in the case illustrated was 2250°K. After 
each interval the filament was cooled suddenly by stopping the filament 
heating current so that the filament would pass rapidly through the range 
of activating temperatures. The temperature was then raised to the 
testing temperature and the emission measured. 

Nature of the thorium film on activated tungsten filaments. The behavior 
of thoriated tungsten filaments suggests that the thorium, when it dif- 
fuses to the surface, forms an adsorbed film which completely covers the 
surface. The maximum activity reached by heating at different activ- 
ating temperatures is the same although the rate of activation may 
vary greatly. It may seem strange at first that thorium can accumulate 
in the surface to an apparently high concentration while the concentration 
in the interior is very small. But this phenomenon is not different from 
that observed in all solutions containing a solute capable of lowering the 
surface energy of the solution. For example ether dissolved in water 
lowers the surface energy and therefore by Gibbs’ rule becomes adsorbed 
in the surface. 
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The properties of thorium are such that we should expect it to lower 
the surface energy of tungsten. The Eétvés-Ramsay-Shields relation 
indicates that the surface energy of any substance at temperatures far 
below the critical temperature should be proportional to T,/V™ where 7, 
is the critical temperature and V is the atomic volume. Critical tempera- 
tures of substances of similar type are proportional to the boiling points 
at ordinary pressures. Thus for different metals the surface energies 
should be proportional to 7,(d/A)** where T; is the boiling point, 
d is the density, and A is the atomic weight. Thorium seems to be about 
as volatile at 2000°K as tungsten is at 3000°K. The density of thorium 
is 11.9 while tungsten is 19.4 and the atomic weights are 232 and 184 
respectively. Thus each of the three factors differs for the two metals in 
such a way as to make the surface energy of thorium less than that of 
tungsten. In fact, from the relation given above, we can calculate that 
the surface energy of thorium should be only 41 per cent of that of tung- 
sten. Thorium dissolved in tungsten should therefore tend to lower the 
surface energy and by Gibbs’ rule should be adsorbed in the surface. 

During the last few years the writer has found evidence of many kinds 
that adsorbed films do not in general exceed a single atom or molecule 
in thickness.* The film of thorium responsible for the high electron emis- 
sion from an activated filament should also not exceed an atom in thick- 
ness. The deactivation curves such as those illustrated in Fig. 1 furnish 
very direct confirmation of this view. The experiments show that the 
shape of the deactivation curve is independent of the length of time 
during which the filament has been allowed to activate. Thus curves 
AC and BD in Fig. 1 are identical although the thorium film was allowed 
to build up for a much longer time in one case than the other. If the film 
had become thicker than a single atom after the long activation at 2100° 
there would have been a delay in the falling off of the emission so that the 
slope of the upper part of BD would have been much less than that of AC. 

There is another indication that the film is never more than one atom 
thick. If the film of thorium should continue to grow in thickness, we 
should not expect the increase in emission to follow any simple law, but 
should find a radical change in the law of increase when the film became 
more than one atom thick. We shall see, however, that the activation 
curves such as OAB in Fig. 1 are in excellent agreement with simple 
mathematical equations over the whole range from a fully deactivated 
condition to that of full activation. 


3A full summary of this evidence with references to the original publications was 
given in Trans. Faraday Soc, 17, 607 (1921) and reprinted in Gen. Elec. Rev. 25, 
445 (1922). 
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With this conception of an adsorbed film of thorium one atom thick, 
we can see that the activity of the surface should depend on the fraction 
of the surface covered by thorium atoms. We shall represent this fraction 
by 6. We now need to determine the relation between @ and the electron 
e mission 17. 

It will be shown in a separate paper that the relation between 7 
and @ may be deduced by means of the electric image theory proposed by 
Schottky.5 With the accelerating fields ordinarily used in measuring 
saturation currents, the electron emission depends upon the number of 
electrons that have velocity components sufficient to enable them to 
reach a certain critical surface, about 5X10-* cm out from the cathode 
surface. At the critical surface the force due to the accelerating field 
is equal to the attractive force (image force) acting on the electron 
because of the positive charge which the electron induces in the metal- 
lically conducting cathode surface. The effect of thorium in increasing 
the emission is not due to a specifically high emission from its surface 
but results from the positive contact potential of the thorium surface 
with respect to the tungsten surface. This positive potential of portions 
of the surface raises the potential at the critical surface and thus increases 
the number of electrons that can pass this critical region. From this 
theory it is found that the logarithm of the emission is a linear function 
of @ so that 6 may be calculated from the equation 

a log 1 —log % 
log 1: —log % 

where % is the saturation current corresponding to pure tungsten (@=0) 

and 7; is that from the fully activated surface (@=1) while z is the current 

when the fraction of the surface covered by thorium has any value 

represented by 8. 

If, for example, the emission from tungsten at a given temperature is 
one unit while that from a fully activated filament at the same tempera- 
ture is 100,000 units, then the emission from a surface half covered with 
thorium (@=0.5) is not the arithmetic mean of 1 and 100,000 but is the 
geometric mean 316. Of course if the thorium on the half covered surface 
should be placed all at one end so that one half the length of the filament 
was completely covered with thorium while the other half had no thorium 
the emission would actually be given by the arithmetic mean 50,000. 
Thus the distribution of the thorium can influence the saturation current. 

The electric image theory of the effect, however, indicates directly 
how uniform the distribution must be in order that Eq. (1) may apply. 





(1) 


’ Schottky, Phys. Zeits. 15, 872 (1914) 
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The thorium influences the emission only by affecting the potential at the 
critical surface which is at a distance from the surface equivalent to about 
200 atomic diameters. Any distribution of thorium which gives a nearly 
uniform potential at the critical surface will be sufficiently uniform to 
cause the emission to vary in accord with Eq. (1). Analysis shows that 
even if the thorium were concentrated into squares like those of a checker- 
board Eq. (1) would apply fairly accurately if the squares were not larger 
than about 100 atoms on a side. With any distribution likely to be 
reached on a plane surface by diffusion or evaporation we should then 
expect Eq. (1) to apply. 

Another relation which follows from the theoretical study of this 
problem is that the value of } in Richardson’s or Dushman’s equation 
for electron emission should also be a linear function of @ so that 

6 = (b—bo)/(b1 — bo) (2) 
where bo and 5; correspond to the values of 6 for pure tungsten and for a 
fully activated surface, respectively. 

The shape of the activation curve. Whenever the thoriated filament 
is heated to a high temperature at which deactivation occurs, all thorium 
evaporates off the surface as fast as it arrives from the interior by dif- 
fusion. In the filament just below the surface the concentration gradient 
should be uniform, when the concentration at the surface is zero. Thor- 
ium should thus diffuse to the surface at a practically constant rate during 
an activation of the filament, except in so far as the supply of thorium 
in the filament becomes exhausted. Since the amount of thorium needed 
for activation is only that in a film one atom thick, it seems probable 
that the available thorium in the filament should suffice for many 
activations. Experiments show in fact that in a series of activations 
alternating with deactivations at rather low temperatures where very 
little thorium is produced by reduction, the rates of activations decrease 
only very slowly in successive runs. During the early part of any single 
activation run we should thus expect the thorium in the surface, and 
therefore 6, to increase at a uniform rate. This means that the logarithm 
of the current should increase in proportion to the time, or the current 
should increase exponentially with the time. This is in general accord 
with the shape of the activation curve shown in Fig. 1, or curve A in Fig. 2. 

Since the logarithm of the current varies linearly with @ it is more 
reasonable to plot an activation curve by using log 7 as ordinate, Curve 
B in Fig. 2 is plotted in this way from the same data as those used for 
Curve A. The current 0.01 is that from a completely deactivated sur- 
face, the dotted line MN thus corresponding to @=0. The current 1000 
is from a fully activated surface, the dotted line S’R’ corresponding to 
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6=1. Distances measured vertically between these lines MN and S’R’ 
are directly proportional to 6. Thus during the initial part of the activa- 
tion run, before the current had increased to a value that can be shown on 
Curve A, the value of 6 had risen from 0.0 to about 0.6 and the current 
had increased about one thousand fold. The fact that the slope of Curve 
B changes so little in this interval is good confirmation of the linear 
relation between log 7 and @. 

If the thorium diffuses to the surface at a uniform rate and remains 
in the surface, @ should increase at a constant rate and log 7 plotted against 
t, should give a curve like MQR’ instead of the observed curve B. It 
is clear therefore that the thorium accumulates at a decreasing rate as @ 
increases. It was at first thought that such an effect might be due to 
diffusion from the highly concentrated surface film back into the mass 
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Fig. 2. Rate of activation of thoriated tungsten filaments. 


of the filament. It would follow from this that the net rate of arrival 
of thorium at the surface should fall to zero as the surface film becomes 
filled (@=1); the concentration gradient near the surface would also 
approach zero, while the concentration itself, just below the thorium 
film, would increase. The thorium which failed to accumulate in the 
surface would thus pile up under the surface. It is an easy matter to 
test whether this is so, for the increased concentration would raise the 
initial rate of activation in a subsequent run. 

Experiments have shown that the thorium does not increase in con- 
centration in this manner after a long activation. Ina typical experiment 
of this kind (Exp. 528, runs M, to P;) the filament was first flashed at 
2800° for 23 sec. to get a high thorium content, then heated to 2500° 
for 100 sec. to allow irregularities in the thorium concentration to 
become smoothed out. The filament was now run at 2055° for 187 
minutes (run M,). The initial rate of activation (corresponding to the 
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straight line MQ in Fig. 2) was dé/dt =0.00060 per sec. while the final 
rate at the end of the run was less than 0.000003. The filament was now 
flashed at 2800° for less than one second to distill the thorium from the 
surface without producing any appreciable fresh supply of thorium. 
On running the filament again at 2055° (run N;) the initial rate of activa- 
tion was 0.00034 showing that the very long previous activation had not 
caused any accumulation of thorium under the surface but had resulted 
in a decrease in the thorium concentration. 

The filament was kept at 2055° (still run N,) for a total time of 277 
minutes, the activation at this time having come to a standstill (rate less 
than 10~° per sec.). In order to make sure that the flashing for 1 sec. 
in the previous run (which had completely deactivated the filament to 
6=0) had not vitiated the results by causing a loss of thorium under the 
surface, the filament was now only partly deactivated by flashing for 3 
sec. at 2500°. This flashing lowered the value of 6 to 0.78. The filament 
was brought again to 2055° (run O,). The rate of activation was now 
dé/dt =27 X10-* per sec. whereas in the previous run (N;) at a similar 
value of @ the rate had been 61 X10~. If instead of continuing the runs 
for such great lengths of time, they were stopped after say 30 or 40 
minutes, the rates of activation fell off in succeeding runs much less than 
in the examples shown here. In no case were any indications found of an 
accumulation of thorium under the surface such as would result from back 
diffusion from the surface film. We may therefore conclude that under 
ordinary experimental conditions the formation of the adsorbed film of 
thorium is an irreversible process in the sense that thorium never goes 
from the surface back into the body of the filament. 

It is clear from these experiments that thorium diffuses continually 
to the surface of a filament even after the surface film becomes saturated. 
Since it does not accumulate in the film it must leave the surface in the 
form of vapor. This is not, however, due to ordinary evaporation of the 
adsorbed thorium film, for it occurs at a much lower temperature. The 
phenomenon we have just been discussing (falling off of d@/dt with time) is 
just as well defined at an activating temperature of 1900° as at 2055° but 
evaporation of a completely formed adsorbed film takes place very slowly 
until temperatures of about 2150° are reached. 

Examination of a large number of logarithmic activation curves (like 
B in Fig. 2) shows that they are accurately represented by an exponential 
formula of the type: 

log (1—0)=—kt or 0=1-e*" (3) 
from which we obtain 
do/dt=k(1—8@). (4) 














a 


7 ees ob a 














368 IRVING LANGMUIR 


Let No be the number of thorium atoms per cm? in a complete mona- 
tomic adsorbed film (@=1), and let G be the concentration gradient just 
under the adsorbed film, expressed in terms of atoms of thorium per cm* 
per cm (atoms Xcm“). If D is the diffusion coefficient (in ordinary c.g.s. 
units, cm? sec.“') of thorium through tungsten at the activating tempera- 
ture, then DG is the rate per unit area at which thorium atoms arrive 
at the surface by diffusion from the interior. If the surface has just 
been deactivated (@=0) no thorium leaves the surface as vapor, but that 
which first reaches the surface accumulates there. The rate of accumula- 
tion of thorium atoms is Nod@/dt and we may therefore equate this to DG 
when @ is very small, that is in the range where the curve B (Fig. 2) 
practically coincides with the straight line MQ. From Eq. (4) we see 
that dé/dt for larger values of @ is proportional to (1—6). We may 
therefore place for all values of @ 

N,d0/dt = DG(1—6). (5) 

By comparison with Eq. (4) we see that DG is directly related to k, the 
rate of activation given by the experiments: 

DG = Nok. (6) 

From Eq. (5) it appears that the rate at which thorium atoms leave 
the adsorbed film in the form of vapor is DG@. It depends not only 
on the condition of the surface which is defined by 6, but upon the rate 
at which the thorium diffuses to the surface. We shall therefore refer to 
this type of vaporization as induced evaporation in order to distinguish 
it from normal evaporation which depends only upon the thorium present 
in the surface film. 

The induced evaporation seems to be caused by the difference in the 
normal rates of evaporation of thorium atoms from tungsten and from 
thorium. Thus suppose that the normal rate of evaporation of separate 
thorium atoms in contact only with tungsten atoms is negligible at the 
activation temperature under consideration, but that the rate of evapora- 
tion of thorium atoms which are held by one or more underlying thorium 
atoms is very high compared to the rate of diffusion of thorium to the 
surface. If a thorium atom diffuses to a point on the surface already 
covered by thorium, it will cause the evaporation of one of the surface 
atoms as soon as the incoming atom reaches a position just under the 
surface atom. The new atom then takes the place of the atom which 
has just evaporated, but has not contributed anything to the thorium 
content of the surface. The only case in which the thorium content 
of the surface does increase is when a thorium atom diffuses to a point 
in the surface which is not already occupied by thorium. Since the chance 
that any given place on the surface is unoccupied is (1—@) we are led 
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directly to the conclusion that the rate of growth of the thorium film is 
proportional to (1—#@) in agreement with the empirical Eq. (4). We have 
thus reached a rational derivation of Eq. (5). 

If in Fig. 2 we measure distances in arbitrary units from any point P’ 
to the dotted line S’R’ and plot the logarithms of these distances as 
ordinates of a moving point P’’ we obtain the line denoted as curve C. 
Curves B and C are conveniently both plotted on the same sheet of semi- 
logarithmic paper. If the experimental observations are in accord with 
Eq. (3) or Eq. (5), curve C will be a straight line whose slope is equal to 
2.303 k, if logarithms to the base 10 are used. 
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Fig. 3. Rate of activation of a thoriated tungsten filament at temperature 7, 
after flashing at 2800°K. (Exp. 527, Y, E:, X:1). Surface of 
filament =0.217 cm?; 7;=1445°K 


Fig. 3 is an illustration of typical experimental data. The saturation 
currents were measured at a testing temperature of 1445°K at intervals 
during the activation runs. The experimentally determined points are 
shown by small circles on the logarithmic plots.A, B and C. Curve A 
is a run at an activating temperature of 2000°K. By trial it was found 
that the limiting value of the current towards which the observed currents 
were approaching asymptotically was 1000 microamperes. A horizontal 
line was drawn corresponding to this current and distances were measured 
vertically in cm from this line to the observed points on Curve A. These 
were plotted logarithmically (by the method described for curve C, 
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Fig. 2) and gave the straight line A’. From this straight line by reversing 
the above process, the curve shown by the full line marked curve A was 
calculated. It will thus be noted that curves A, B, and C in Fig. 3 are 
not merely smooth curves drawn to fit the observations but they are true 
exponential curves calculated from equations like Eq. (3) by the choice 
of suitable values of k. It is seen from Fig. 3 that the observed points 
fall on the calculated curve with high accuracy except for very low values, 
in the neighborhood of log 7;=0, which corresponds to currents of about 
one microampere. Since the currents in this particular set of experiments 
were measured with a portable microammeter the agreement is still 
within the error of measurement. 

Other experiments made at higher testing temperatures or with more 
sensitive instruments have shown that the agreement with the ex- 
ponential equation is, in general, equally good down to the lowest states 
of activity (where 6=0). These results not only confirm Eq. (5), but 
check the accuracy of Eq. (1) and indicate that the whole activation 
curve follows a single simple law. It is this evidence that was referred 
to previously in connection with the proof that the thorium film is 
monatomic. 

Evaporation of thorium and deactivation curves. Although the rates of 
activation could be made to vary widely, even at a given activating tem- 
perature, by changing the concentration gradient G, practically the same 
limiting value of 7,,=1000, at 7;=1445°K, was found repeatedly in 
separate runs at activating temperatures ranging from 1900° to 2100°. 
But if G had an unusually low value, as for example after heating the 
filament for many hours at 2300°, the limiting currents found with 
activating temperatures of 2000° or more were distinctly lower. This 
effect is due to normal evaporation of thorium from the surface at a rate 
comparable with the rate of supply of thorium by diffusion from the 
interior. 

Curve C in Fig. 3 at an activation temperature of 2150° illustrates the 
effect of this evaporation. The rate of activation is much greater than 
in A, for the rate of diffusion is increased by the higher temperature. 
But the limiting current 7, has fallen to 253 microamp., which cor- 
responds to 6., =0.88. 

If we are to take into account normal evaporation of thorium, Eq. (5) 
must be modified by subtracting from the second member a term E 
which we may define as the rate of evaporation expressed in atoms per 
second per cm? of filament surface. The quantity E will depend on the 
temperature and on 6. The most natural supposition, perhaps, is to 
assume that E is proportional to @, viz., the evaporation is proportional 
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to the amount of thorium on the surface. But we shall see that the experi- 
ments indicate that this is not the case. When the surface is very 
largely covered with thorium it is in fact to be expected that the thorium 
atoms will tend to hold on to each other and so make the rate of evapora- 
tion increase less rapidly than in proportion to 6. Since it is convenient 
to have an expression for the variation of E with @ let us assume pro- 
visionally that EZ is a 'inear function of 6. This will always be permis- 
sible within a restricted range of values of @6. We may thus place 
E = vot v0. (7 ) 
Subtracting this from the second member of Eq. (5) gives 
N,od6/dt = DG(1 —0) — (vo +8). (8) 
When t=, d6/dt=0 and 6=8,,, so the equation becomes 
vot+v0,. = DG(1—8., ). (9) 
Subtracting (9) from (8) and integrating 


log (= = ) = —hkt (10) 
6.. —% 


where @ is the value of @ corresponding to ¢=0 and & is a constant 
k=(DG+»)/No. (11) 
The limiting value of @ is @,. given by 
6.. = (DG —%)/(DG+yv). (12) 

The experimental data furnish us with values of k and @,, and from 

these we wish to find DG, » and v. From the above equations we find 
DG = Nok —v = Nok0.. +0 (13) 
votv = Nok(1—8..). (14) 

When @=0 and @,,=1, Eq. (10) reduces to Eq. (3). When the rate 
of evaporation » is not negligible compared to the rate (DG) at which 
thorium diffuses to the surface, we see from Eq. (12) that @,, is less than 
unity, so that the filament does not become fully activated. The linear 
relation between log (1—@) and ¢ given by Eq. (3), which agreed with 
experiments where no evaporation occurred, must now be replaced ac- 
cording to Eq. (10) by a linear relation between log (6,, —@) and ¢. 

Thus in Fig. 3 the lines B’ and C’ were obtained by plotting as ordi- 
nates the logarithms of the vertical distances from points on the curves 
B or C to the horizontal asymptotes corresponding to the values of 7,, 
for the curve in question. It is seen that the agreement between the 
observed points on curves B and C is as good as it was in the case of curve 
A. It should be kept in mind that the curves A, B and C were calcu- 
lated by means of Eq. (10), choosing of course in each case values of 6. 
(or 7.) and k which gave best agreement with the experiments. Activa- 
tion curves taken under the most varied conditions, even at such high 
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activating temperatures that 6,, did not exceed 0.1, always gave reason- 
ably good agreement with straight lines when plotted in the manner of 
curves A’, B’ and C’. 

Run X;, in Fig. 4, is an example of a deactivation curve obtained with 
the filament at a deactivating temperature Tz of 2300°. Preceding this 
run the filament had been flashed 30 seconds at 2800°, then run 3 hours 
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Fig. 4. Deactivation of a thoriated tungsten filament at 7,=2300°K. 
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at 2250° and 3 hours at 2100°. This had caused about 4/5 of the thorium 
produced at the high temperature to diffuse out of the filament and had 
thus lowered G to about 1/5 of its value after flashing. In run X,, the 
activity falls rapidly at first, reaches a minimum and slowly rises to a 
steady value. This is characteristic of all deactivation curves in which 
G has previously been brought to a low value. 

The reason for the minimum is that the activity depends upon two 
factors that act in opposite directions; the activity tends to decrease 
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because of evaporation but the increased production of thorium in the 
filament due to the higher temperature increases G and tends to raise the 
activity. Since the latter effect occurs more slowly than the former it is 
readily seen that the activity should have a minimum. 

A second type of deactivation curve, illustrated by run Sj, Fig. 4, is 
obtained if the filament is flashed at a high temperature so as to give a 
high value of G and is then immediately run at a lower temperature. 
In the case of run S, the filament was flashed for 30 seconds at 2800° 
and then run at 2300° for 106 minutes, the activity being measured 
at intervals at a testing temperature of 1630°. In this case the 
activity rises rapidly to a maximum and then falls slowly to the 
same limiting value (D in Fig. 4) as in the first type of activation 
curve. There is no minimum activity in this case. The reasons 
for this shape of curve are similar to those of the first type, but the 
displacements from the stationary state are in opposite directions. At 
the beginning of the run there is no thorium on the surface so that 
evaporation cannot occur, but the thorium arrives at the surface at a 
rapid rate because of the high value of G. The activity thus increases. 
In the meantime, however, the thorium is diffusing out of the filament 
faster than it is being produced so that the value of G steadily but slowly 
decreases causing finally a decrease in the activity. 

In the latter part of both types of curve of Fig. 4 the thorium evap- 
orates from the surface at nearly the same rate as it arrives from the 
interior so that in a sense a “‘steady” state may be said to have been 
reached although this state gradually changes as G increases or decreases. 
In the first parts of the curves, however, the rates of evaporation and 
diffusion are not even approximately balanced, the conditions approach- 
ing those assumed in the derivation of Eqs. (10) and (11). The steady 
state corresponds to the case where 6=86,, as given by Eq. (12). The 
dotted curve AE in Fig. 4 is drawn as a continuation of EG and the whole 
curve AEG corresponds to values of @,,.. The sharply rising initial part 
of run S;, marked CE, is thus an activation curve like those illustrated in 
Fig. 3 except that the value to which the curve approaches (correspor!ing 
to @.,) is not constant but decreases according to curve AEG. Similarly, 
CBD represents the limiting current values corresponding to 6,, for run X,. 

From the curves of Fig. 4 we see that the times required for the activity 
to go 90 per cent of the way to the “steady” value was 4.5 minutes for 
run S,; and 19 minutes for run X;, these being in a ratio of 1 : 4.2. From 
Eq. (11) this must also have been the ratio of the corresponding values of 
(DG++y) so that DG in run X, must have been at least 4.2 times as great 
as in run S}. 
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We shall return to a quantitative consideration of activation and deac- 
tivation curves after the following section. 

Arrangement of tungsten and thorium atoms in surface films. Tungsten 
crystallizes with a body-centered cubic lattice of atoms, the side of each 
cube being 3.150 A. The shortest distance between atoms (from the 
center atom to an atom in a corner) is /(?) X3.150 or 2.728 A. Ina 
plane (100) through any atom parallel to a face of the elementary cube 
the number of atoms per cm? is 1.008 X10". In a plane passing through 
the center atom and any of the edges of the cube (i. e. 110 plane) there 
are 1.425 X10" atoms per cm’. 

By special heat treatment a square tungsten rod (8X8 X100mm) 
was brought into the form of a single crystal. As it had been held for a 
long time near its melting point, the surface was distinctly etched 
by evaporation so that it reflected light in several directions which 
were not related to the original faces of the rod. By mounting a fragment 
of this rod on a cork and holding the specimen in a dark room so that 
the etched crystal faces reflected light from a single lamp placed almost 
behind the observer’s head, it was possible to place pins in the cork 
indicating the directions of the normals to the crystal faces. Inspection 
of these pins proved that the faces developed by etching are those of the 
rhombic dodecahedron (110 planes), all these 12 faces and no others 
being observed. These are the faces which contain the greatest possible 
number of atoms. Other single crystals of tungsten were chemically 
etched by an alkaline ferricyanide solution but only the same faces were 
found.’ Tungsten surfaces which have been heated to such high tem- 
perature that some evaporation has occurred, thus consist exclusively 
of dodecahedral faces. Simple considerations show that the extent of 
the total developed surface does not increase as the size of the faces 
increases. 

Consider, for example, a single crystal of tungsten having the external 
form of a sphere, and let this be etched until the original surface is 
gone and has been replaced by etched dodecahedral faces. Imagine 
that the sphere be viewed from a point far out along one of the cubic 


® Hull, Phys. Rev. 17, 576 (1921) 

7 The cleavage surfaces of the single crystal rod were found to be cubic faces (100 
planes). Chromium which belongs to the same group in the periodic table, and crystal- 
lizes with the body centered cubic lattice, showed pronounced cubic cleavage (100 planes) 
with a less well developed secondary cleavage along the dodecahedral faces (110 planes). 
After placing a fragment in dilute hydrochloric acid for a few hours the cubic faces 
disappeared and the dodecahedral faces were all brought out distinctly. Large crystals 
of iron, (transformer iron containing a few per cent of silicon) which also showed the 
body centered lattice, gave on etching in acid only the cubic faces (100 planes). 
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axes, say the X-axis. From this position only faces parallel to four out 
of the 12 dodecahedral faces can be seen and each of these faces is inclined 
at 45° to the line of sight. Every part of the crystal that can be seen 
must consist of just such faces. Thus if 7 is the radius of the sphere, the 
apparent or projected area is mr’? but the actual area seen is 4/2: rr’. 

As this represents only 4 of the 12 faces, the total surface of the etched 
sphere is 34/2-xr*? as compared with 47r? for the original sphere. The 
effect of etching has thus been to increase the surface 3/(2./2) = 1.0607 
fold. The depth of etching has no influence on this ratio as long as the 
projected area of the sphere is not materially changed. The total surface 
of a tungsten filament consisting of small crystals oriented at random 
is thus 6 per cent greater than its apparent surface. 

Thorium crystallizes in a face-centered cubic lattice’, the side of the 
cube being 5.04 A. Since this is a close-packed type of lattice, the 
thorium atoms behave much like spheres of diameter 3.564 A. The 
number of such spheres that can be packed per cm? into a plane lattice 
is 0.90910" which is only 64 per cent as many as the atoms in the 
110 planes in tungsten. We shall see from the data on evaporation of 
thorium, that forces between adsorbed thorium atoms and the under- 
lying tungsten atoms are very large while those between thorium atoms 
in adjacent positions on the surface are relatively small. It is therefore 
highly probable that there is a stoichiometric relation between the 
number of thorium atoms in a fully active surface and the number of 
underlying tungsten atoms. A one-to-one relation would involve ex- 
cessive crowding of the thorium atoms and would be incompatible with 
small forces between adjacent atoms. We shall therefore assume that 
there are half as many thorium atoms as underlying tungsten atoms in a 
fully active surface. Thus the maximum number of thorium atoms in 
a plane surface is 0.713610" per cm*. In the case of the surface of 
an ordinary tungsten filament we must apply the correction due to 
etching (1.0607) so that the quantity that we have called No in Eqs. (5) to 
(14) should have the value 

9 =0.756 X10” atoms per cm’. (15) 

Rate of activation as a function of temperature. At temperatures so 
low that evaporation of thorium is negligible, the’ rate of activation, as 
measured by &, is proportional to the product DG according to Eq. (6). 
By proper preliminary heat treatment G can be brought to the same 
value in each of a series of runs carried out at different temperatures. 
The rate of activation k is then a relative measure of the diffusion 
coefficient D. 


§ Hull, Phys. Rev. 18, 88 (1921) 
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Table I gives the results of such a series of runs (in the latter part 
of Exp. 527). Immediately preceding each run the filament was flashed 
at 2800°K for 5 seconds to bring G to a definite and constant value. The 
filament temperature was then changed to Tq (Col. 2), this being inter- 
rupted at intervals ranging from 5 to 20 minutes for the purpose of 
measuring the activity at the testing temperature (7;) of 1445°K. The 
saturation currents at 7; were plotted on semi-logarithmic paper as 
illustrated in Fig. 3. The limiting current 7, corresponding to t= is 
given in microamperes in Col. 3. The filament surface from which 
these currents were obtained was 0.217 cm*. By plotting logit, — 
logit as ordinate (on the same sheet of semi-logarithmic paper), against 
time as abscissa, straight lines were obtained as in curves A’, B’, and 
C’ in Fig. 3. The slopes of these lines, after multiplying by 2.303 to 
convert to natural logarithms and dividing by 60 to convert from rates 
per minute to rates per second, give directly the values of k defined by 
Eq. (10). These are tabulated in Col. 4. Col. 5 gives 0,, calculated by 
Eq. (1) from i,, of Col. 3 using the values 1;=1300 and i=0.013 
microamperes. Col. 6 gives values of 0, calculated by means of Eq. 
(14) from data on » and » and will be discussed in a later section. 

TABLE I 


Rates of activation at various temperatures 
Experiment 527 Filament surface 0.217 cm? 





(1) (2) (3) (4) (5) (6) 
Run i k 


a lo 


Oo On 
w-amps. at (10-4 per sec.) (obs.) (calc.) 
Ti = 1445°K 





2050°K 790 5.31 0.957 0.958 
2050 738 6.36 0.951 0.964 
2150 253 15.41 0.856 0.856 
2105 477 9.29 0.913 0.913 
2050 707 5.11 0.948 0.957 
2050 805 5.15 0.958 0.957 
1950 1300 1.70 1.000 0.990 
2000 1000 2.59 0.977 0.976 





The diffusion of one metal through another must involve interchanges 
in position between neighboring atoms in the space-lattice. Only atoms 
containing a certain critical amount of kinetic (or potential?) energy 
are capable of making such interchanges. The frequency of these changes 
of position and therefore the diffusion coefficient should thus be pro- 
portional to e@p/R7, where Qp may be taken as the “heat of diffusion” 
or the critical energy (per gram atom) involved in atomic interchanges.® 
From this relation it appears that the natural logarithm of the diffusion 


*S. Dushman and I. Langmuir, Phys. Rev. 20, 113 (1922) 
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coefficient is a linear function of the reciprocal of the temperature, the 
slope of this line being equal to Qp/R. 

To test this conclusion the values of logik from Table I have been 
plotted, in Fig. 5, against 1/T,. Since 4 runs were made at T7,=2050°, 
the average was taken of the values of k for these runs, viz. k = 0.000548. 
The points lie close to a straight line. Taking R=1.987 calories this 
slope gives Qp = 94000 calories per gram atom. 

Absolute value of the diffusion coefficient. To determine the actual 
value of the diffusion coefficient we need to measure some property 
involving the total thorium content of the filament. Heating the filament 
momentarily to a very high temperature increases the concentration of 
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Fig. 5. Rate of activation of thoriated filaments as function of temperature. 


thorium. By heating the filament for a long time at a much lower 
temperature this supply of thorium gradually diffuses out so that the 
concentration gradient at the surface decreases in proportion. By 
measuring the product DG from time to time during this exhaustion 
of the thorium supply, it is possible to determine the absolute value of D. 

The diffusion problem involved is analogous to the thermal problem 
of the cooling of a long cylinder suddenly plunged into a perfectly 
conducting fluid maintained at a fixed temperature. 

This problem has been treated by Carslaw’® who gives the general 
equations for the temperature distribution in an infinite cylinder resulting 
from any initial distribution. For the particular cases required in the 


1° Carslaw, Fourier’s Series and Integrals, Macmillan and Co., 1906, p. 316 
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present problem Mr. Harold Mott-Smith has derived the following 
equations. 

Consider an infinitely long circular cylinder in which the concentration 
(of thorium) is a function f(7) involving only the distance r from the 
axis of the cylinder. At a certain instant, which we shall choose as 
origin of time, (‘=0), the concentration over the whole surface of the 
cylinder is brought to zero and kept at zero. This is actually accom- 
plished by raising the temperature of the filament suddenly at the time 
0 to such temperature that all the substance that diffuses to the surface 
evaporates immediately. Then after the time ¢, the concentration 
gradient G at the outer surface of the cylinder is expressible as an infinite 
series of exponential terms: 

G= Be? (16) 
where 
P= u,2Dt/R? (17) 
Here R is the radius of the cylinder and D is the diffusion coefficient. 
The quantities u, are the successive roots of the equation Jo(x) =0, where 
Jo(x) is the zeroth order Bessel function of x. In using the foregoing 
equations it is merely necessary to substitute the following numerical 
values of these constants :" 
ur= 5.783 Ms’= 74.89 
po? = 30.47 ba? = 139.04 (18) 

The coefficients B depend upon the initial distribution of the diffusing 
substance, at the time ¢=0. In case the concentration is initially uniform 
throughout the cylinder and equal to C, we find 

B,=B,.=B,;= .... =2C,/R (19) 

If the diffusing substance is being produced at a uniform rate through- 
out the whole body of the cylinder while the concentration at the surface 
is kept zero, the concentration finally assumes a parabolic distribution 
given by 

C=C,|1—(r/R)?] (20) 
where C>» is the concentration at the axis of the cylinder and C is the 
concentration at any point whose distance from the axis is r. The 
concentration gradient G at the surface is then 

G=2C,/R (21) 

If the production of diffusing substance is suddenly stopped after 
this parabolic distribution has been reached, as for example by lowering 
the temperature, the subsequent changes in G are still given by Eq. (16) 
but the values of the coefficients B now become 

B,=8C)/(un?R) (22) 


11 Smithsonian Physical Tables, 1920, p. 68 
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From Eq. (19) the initial value of G is infinite for the case of initial 
uniform distribution. With the parabolic distribution, we see from Eqs. 
(16) and (22) that the initial value of G is 8C,[2(1/y,2)]/R which is the 
same as given by Eq. (21). 

From the rapid increase in the exponents (18) in the series (16) it is 
evident that no matter what the initial distribution may be, the higher 
terms of the series soon become negligible compared to the first term. 
Curve II in Fig. 6 is a logarithmic plot of G as a function of time for the 











“00! 2 003 004 005. 008 01 02 03 04.050 06 1 2 ~. : 2 


Fig. 6. Transient effects in diffusion in filaments 
Curve I. initial uniform distribution in semi-infinite solid with plane face. 
Curve II. Initial uniform distribution in cylinder. 
Curve III. Initial parabolic distribution in cylinder. 
Curve IV. First term in series, Eq. (23). 


case of initially uniform distribution corresponding to Eqs. (19) and 
(16). As ordinate, RG/C, is plotted on the logarithmic scale, while 
Di/R? is the abscissa on the same scale. This method of plotting has 
the advantage that any factorial change in R or D or the units in which 
these or ¢ are expressed, does not change the shape or scale of the curve, 
but is equivalent to a displacement of the origin of co-ordinates. Curve 
IV corresponds to the first term of the series (16) with the value of By; 
from Eq. (19): G=(2C,/R)e™ (23) 
where by Eqs. (17) and (18) 
a=5.783D/R? (24) 
It is seen from Fig. 6 that even when the initial concentration is 
uniform, G is given with high accuracy by Eq. (23) for all values of ¢ 
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which make Dt/R? greater than about .16. Let us consider on the other 
hand the case where the diffusion has proceeded for only a short time so 
that the concentration in the cylinder differs appreciably from the original 
uniform concentration only in layers whose depth below the surface 
is small compared to the radius. Then the conditions are essentially 
similar to those occurring when the plane surface of a semi-infinite 
solid is suddenly changed in temperature, a problem treated in most 
books dealing with heat conduction.” Thus for the corresponding dif- 
fusion problem, when the initial concentration is uniform and the con- 
centration at the surface is brought to zero, we find 
G=C,/VxDt (25) 
To compare this with the functions plotted in Fig. 6 it is conveniently 
transformed into 


RG/C,,=0.564-V R2/(Dt) (26) 
which corresponds to the straight line I. For very small values of 
Di/R? Curve II approaches I but rather slowly. 

If the initial distribution is of the parabolic type given by Eq. (20), the 
concentration gradient G is also expressible, for sufficiently large Dt/R?, 
by the first term of Eq. (16) which now by Eggs. (17), (18), (22) and (24) 
takes the form | 

G=1.3833(C)>/R)e™. (27) 

By comparing this with Eq. (23) we see that both these functions 

can be expressed by Curve IV in Fig. 6 if we take 
Cp =1.4458C,. (28) 

In other words the distribution finally becomes the same whether 
we have initially a uniform concentration or a parabolic distribution 
with the concentration at the axis 1.446 times as great. Curve III in 
Fig. 6 gives the values of RG/C, for this case calculated by means of 
the complete series of Eq. (16) using the coefficients given by Eq. (22). 
Eq. (27) thus gives a good approximation whenever Dt/R? exceeds 0.1. 

The total amount g of diffusing substance in the cylinder per unit of 
length is tR°C, when the concentration is uniform. For the parabolic 
distribution we find from Eq. (20) by integration 

dp =3aR°C>. (29) 

When the distribution has become such that G falls off in accord with 
Eq. (23) or with Curve IV in Fig. 6, we may find g at any time ¢ from 
the value of G in Eq. (23) by calculating the total amount of substance 
that diffuses out to the surface of the cylinder between ¢=# and t=. 
This gives: 

q= (4rR°C,/ ws) % =2.1729R°Cye (30) 


12 See for example Ingersoll and Zobel, Theory of Heat Conduction, p. 77 
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The distribution within the cylinder under these conditions (i. e., 

when G is expressible as the first term of the series of Eq. 16) is given by 

C=0.8018RGJo( u:r/R) (31) 

where C is the concentration at any point whose distance from the axis 
is r. The concentration at the axis is thus equal to 0.8018 RG. 

The absolute value of the diffusion coefficient is most easily found 
from experiments carried out under conditions such that G varies with 
time in accord with Eqs. (23) or (27). The experiments on rates of 
activation of the filament in successive stages thus furnish data for the 
evaluation of a and from this by Eq. (24) the value of D is calculated. 

Several sets of experiments were made to determine D at various 
temperatures. Most of these were carried out before the complete 
theory of the diffusion had been worked out, so that the experimental 
conditions have not always been chosen as well as would otherwise 
have been possible. 

Table II contains data for the calculation of the diffusion coefficient 
at 2500°K, obtained from Exp. 528. Preceding each of these runs 
I, to M, the filament was flashed at 2800°K for 23 seconds in order to 
bring the thorium content to a definite value. The filament was then 
maintained at 2500°K for the time given in Col. 2 with the object of 
causing the concentration gradient G to decrease according to Eq. (16) 
After the heating at 2500°, the rate of activation k(Col. 4) was measured 
by running the filament at 7,=2055°K, interrupting this run at frequent 
intervals for testing the activity at 7;=1450°K in the usual way. 


TABLE II 
Rates of activation at 2055°K following deactivation at 2500°K for various times 
Experiment 528 Filament surface 0.217 cm? 








(1) (2) (3) (4) (5) (6) (7) 
Run Time to k 0 0x k—kr 
at 2500°K at 1450°K at 2055° fromio (Eq. 42) 
(sec.) (u-amp.) (per sec.) (obs. ) (calc.) (xX 10-*) 





10 760 11.49 x 10-4 0.976 0.981 7.85 
600 548 4.31 0.947 0.949 0.67 
10 780 9.58 0.978 0.977 5.94 
3660 505 3.64 0.941 0.940 0.00 
100 720 6.22 0.970 0.965 2.58 











The values of @,, in Col. 5 were calculated by Eq. (1) from i,, in Col. 3 
taking ip =0.01 and 7;=1000 microamp. The calculated @,, in Col. 6 will 
be discussed later. 


Examination of the data in Col. 4 shows that k decreases rapidly at 
first and then approaches a limiting value of 3.64 10-* which we shall 
represent by k,. In this stationary state thorium is being reduced from the 
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oxide at the same rate as it is lost by diffusion and evaporation. The 
values of k are proportional to G, so that the equations that we have 
derived for the variations of G should be directly applicable to our ex- 
perimental data. Because of the constant production of thorium, how- 
ever, we use for this purpose not k but k—k, as given in Col. 7 of Table IT. 
Thus after sufficient time has elapsed for the effect of the higher terms 
in Eq. (16) to die out, the decay. of the quantity k—k, should follow the 
exponential law given by Eq. (23). The following method allows us to 
compare the experimental data directly with the complete Eq. (16) 
as represented by the curves in Fig. 6. 

It is evident that a plot of log(k—&,) as a function of log ¢ should have 
the same shape and scale as the corresponding curve in Fig. 6 and should 
differ only by a parallel displacement. The points marked J;, Ki, M., 
and J; are placed in Fig. 6 by plotting 10'(k—k,) from Col. 7 of Table II 
against 10-°-¢ where ¢ is the time in seconds fram Col. 2. By trans- 
ferring these points to tracing paper and sliding the paper so as to keep 
its lower edge horizontal, it is possible to bring any two of the points 
such as J; and M, to lie on any desired curve. For example by displacing 
the tracing paper parallel to the line J,A the points J; and M, are brought 
on to the curve III, which corresponds to the case of parabolic distribu- 
tion. The points J; and K,, however, do not fall on this curve but lie 
between II and III. By displacing the tracing paper by the distance 
AB the points J; and M, are brought on to the curve II, corresponding 
to uniform initial distribution. The points J; and K, then lie too low. 
It appears from this that the distribution of thorium is intermediate 
between the two cases considered. It is more likely, however, that the 
true distribution resulting from the production of thorium by heating 
the filament to high temperature approaches closely to the parabolic 
case if sufficient time is allowed for a stationary state to be reached. 
It is rather doubtful if the 23 seconds allowed in these experiments was 
long enough. Failure to reach the stationary state would cause deviations 
of the type observed. 

By having placed the observed points on the theoretical curve in this 
manner, we determine the absolute value of Dt/R?, and from this since 
t and R are known, we find D. Thus if we assume initial parabolic 
distribution we find Dt/R? for the point A to be 0.27. For this point 
t=600 sec., and therefore since R=0.00389 cm we find D=6.8 X 10°° cm? 
per second for 2500°K. Similarly if we assume initially uniform dis- 
tribution, Dt/R? for point B is found to be 0.18 which gives D =4.3 X10-° 
cm?/sec. Of course the same values for D are found from the other 
pairs of points. 
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The diffusion coefficient at 2055°K was determined by a series of 
prolonged activation runs,. N,, O,, and P; carried out at this temperature. 
Just preceding run JN, the filament was flashed for a couple of seconds 
at 2800°K. In run JN, the filament was maintained at 2055° for 277 
minutes, the activity being tested at frequent intervals in the usual way. 
The filament was then flashed at 2500° for 3 seconds, a time so short 
that the filament was only partly deactivated and G was not altered. 
Run O, then followed, the filament being kept at 2055° for 110 minutes. 
After flashing for 4 seconds at 2500° an activation curve was again made 
(run P,) at 2055°. Selecting from these runs points corresponding to 
a definite value of @ we obtain the following data” 


Run t k | k—k, 
N,; 30 min. 3.59X10  3.0610~ 
0; 305 1.62 1.09 
P, 475 1.25 0.72 


Plotting these on tracing paper and comparing with the Curve III of 
Fig. 6 the most probable value of the diffusion coefficient at 2055° is 
found to be 1.1 10-!° cm? sec."!. 

The deactivation curve of run S; shown in Fig. 4 furnishes excellent 
data for the calculation of D at 2300°. We have seen that the portion 
EG of the curve CEG results from the gradual decrease in G as the supply 
of thorium becomes exhausted. In the region represented by CE the 
rate of change of activity is so great that the term Nod@/dt in Eq. (8) is 
at least comparable in magnitude with the terms in the second member 
of the equation. In the region EG, however, d0/dt has become so small 
that the corresponding term may be neglected in comparison with those 
in the second member. Eq. (8) thus becomes 

DG = (»+v0)/(1—6@) (32) 

The quantity G given by this equation cannot be used directly to 
determine the diffusion coefficient according to Eq. (16), for the constant 
production of thorium at 2300° must be allowed for. Let G, be the limit- 
ing value of G which would be reached when the accumulated surplus 
of thorium has been exhausted, the rate of loss of thorium by diffusion 
being then equal to the rate of production. The value of G, can be 
expressed in terms of 6,, the corresponding value of 6, by an equation 
similar to Eq. (32). We now substitute in the diffusion equations the 
quantity (G—G,) just as we previously used (k—k,). Eq. (32) gives 

D 1 1 


(G—G,) =—— — 
v+vo 1-6 1-06 


13 The value of k, was taken to be 0.53 X 10-4 as calculated from values of G, which 
will be discussed later. 


(33) 
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The quantities in the second member are both easily found from the 
experiments. Taking the logarithm of the second member and plotting 
against the logarithm of the time as in Fig. 6, the absolute value of D is 
found. The data of run S; in Fig. 4, are in complete agreement with 
the curves of Fig. 6, in the range beyond ¢=40 minutes corresponding 
to the region where the curves II, III and IV coincide. For smaller 
values of #, the experimental curve lies somewhat above Curve II. 
This should not be taken as evidence, however, that there is a uniform 
rather than a parabolic distribution of thorium in the filament for it is 
probable that the deviations from the Curve III are due to variations 
of the rate of evaporation as a function of 6, from the linear relation 
assumed in Eq. (7) upon which Eq. (33) is based. The assumption of 
this linear relation is always justifiable over a narrow range of values 
of #, but becomes open to question when applied to the wide range 
covered by the data of run S$}. 

If we therefore confine our attention to the portion of run S; taken 
after the first 40 minutes we find by our comparison with Fig. 6 the 
value D=1.1210-* cm? sec. Since in this region the effects of the 
uncertain initial distribution of thorium are eliminated, so that only the 
first term in the series of Eq. (16) need be considered, a simpler method 
may be employed. Plotting the logio of the second member of Eq. (33) 
against the time ¢ (on semi-logarithmic paper), a straight line is obtained. 
The slope of this line, multipled by 2.303, gives directly the value of a 
defined by Eq. (24). From this the value of D is found, and is of course 
the same as that obtained by the curves of Fig. 6. 

Run 7, was made to determine the diffusion coefficient at 2400°. 
The filament was first flashed at 2800° for 30 seconds and then run at 
2400° for deactivation, the activity being measured at a testing tempera- 
ture of 1800°. The run was interrupted three times for the purpose of 
determining the value DG/No by means of very short activation runs of 
5 to 10 minutes at 2100°. The values of @ and d@/dt were found in 
these runs from which, by Eq. (8), DG/No could be calculated, the term 
(vo+v0) being applied as a small correction determined from data which 
will be discussed later. The results obtained from run 7) are 








t DG/No D(G—Gr)/No 
(minutes at 
2400°K) (at 2100°K) (at 2100°K) 


4.7 7.67 X 10-4 sec.-! 3.86 X 10-4 sec.-! 
23.1 4.65 0.84 
45.6 3.95 0.13 
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All three of the observations recorded in this table, when compared 
with the curves of Fig. 6, are found to lie in the region where the Curves 
II, III and IV coincide. It was therefore assumed that the exponential 
type of equation such as (23) is applicable. The value of G, was found’ 
by trial so as to make the three points corresponding to the data in the 
third column give a straight line when plotted on semi-logarithmic 
paper against the time ¢#. The slope of this line gave a=0.00140 
sec.“!, from which Eq. (24) gives D=3.57 X10-® cm? sec."! at 2400°K. 

Table III contains a summary of the absolute values of the diffusion 
coefficient obtained in the way just described. The logarithms of 
these values of D are plotted against the reciprocals of the temperatures 
in Fig. 7. The straight line marked D has been drawn to have the 
same slope as the line in Fig. 5, which corresponds to a heat of diffusion 
of 94000 calories per g-atom. The agreement with the observed points 
is very satisfactory, especially when it is considered that several different 
methods have been used in determining the relative and the absolute 
values of D. 


TABLE III 


Experimentally determined values of the diffusion coefficient of thorium through tungsten 








Run Temperature Diffusion coefficient D 





MN, to P, 2055°K 
Si 2300 
Ti 2400 

I, to M, 2500 











The straight line in Fig. 7 gives the following equation for the diffusion 

coefficient at any temperature (in cm’* sec.~'): 
logio D=0.044 —20540/T (34) 

Table IV gives values of D calculated by this equation. 

Rate of production of metallic thorium. The data furnished by low 
temperature activation curves enable us by Eq. (6) to calculate DG and 
now that we know the values of D from Table IV we can obtain the 
value of G. We have seen from Table II that after the filament had 
been heated a long time at 2500°, k, (at 2055°) approached a limiting 
value of 0.000364. From the value D =1.12X10-? and by Eqs. (6) and 
(15) we thus find G,=2.46X10% atoms of thorium per cm‘ at 2500°. 
In a similar way from data of about 20 other sets of activation runs, 
values of G, for temperatures ranging from 2055° to 3000° have been 








386 IRVING LANGMUIR 


obtained. The logarithms of these values plotted against the reciprocals 
of the temperatures give a straight line whose equation is 
logio G, = 25.22 —9620/T. (35) 
The slope of this line multiplied by 4.575 gives Q=44000 calories 
per g-atom. Let , be the rate of production of thorium per unit volume 
throughout the filament. When the steady state has been reached in 
which the concentration gradient at the surface is G,, the rate of pro- 
duction of thorium per unit length is 2tRDC, and thus 


pr =2DG,/R (36) 
TABLE IV 


Diffusion coefficient of thorium in tungsten and the rate of evaporation of thorium adsorbed 
on tungsten 








E T D E 
(atoms per (cm? sec.-') (atoms per 
cm? sec.) cm? sec.) 


(cm? sec.-') 





1400°K 


6xX10-" 0.445 


2200°K 5.10107! 1.61 X10" 


1500 
1600 
1700 
1800 
1900 


4xX10-" 


15 x 1078 
9x10" 
71X10" 


58.5 
4.18 X 10° 
1.81 X 10° 
5.15 x 106 

03 X 108 


2300 
2400 
2500 
2600 
2800 


1.30 x 10-° 
3.06 X 10-° 
6.73 X10°° 
1.39 x 10-8 
5.11 10-8 


1.22 X10" 
7.80 X 10"? 
4.31 10% 
2.08 X 10" 
3.48 X 10" 


2.3 

2.2 

1.61 X 10-8 
9. 

4.2 

1.7 

5.9 

1 


2000 
2100 


1.58 X 10-7 3.99 x 10% 








1 
4X10" 1.53 X 10° 3000 
.83 X 10719 1.75 x 10"° 








The temperature coefficient of p, is thus the sum of those of D and G,. 
Substituting into Eq. (36) the values of D from Eq. (34), G, from Eq. 
(35) and R=0.00389 cm, we obtain 

logio pr = 27.98 — 30160/T, (37) 
pf, being expressed in atoms of thorium per sec. cm*. The heat of the 
reaction involved in the production of the thorium is thus Q, = — 138000 
calories per g-atom, the negative sign indicating that the reaction 
involves the absorption of energy. 

Although in Exps. 527 and 528 the filaments were heated for several 
hours per day for more than 3 months there was no indication of a 
falling off in the rate of production of thorium. Because of the induced 
evaporation, all the thorium produced at the rate p, leaves the filament 
even at relatively low temperatures. The data used in calculating p, 
in Eq. (37) were obtained with a filament containing 1.0 per cent of 
thorium oxide, corresponding to 4.4X10*° atoms of thorium per cm*. 
If the thorium were produced at a constant rate p, the life of the thoria 
would be 4.4 X 10?°/p, seconds or 1.22 X10!"/p, hours. In this way we 
calculate that the thoria supply should last 94000 hours at 1900°. Col. 5 
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of Table VIII contains data on the life at different temperatures. Con- 
sidering that the rate of production of thorium must decrease as the 
thoria content falls off, the ‘‘life” calculated in this way probably repre- 
sents the time required for the thoria to fall to one eth of its initial value. 
For different filaments, even containing the same amount of thoria, 
p, varies considerably depending on the mechanical state of subdivision 
of the thoria in the filament and on other factors. The presence of even 
traces of carbon tends to reduce the thoria to thorium and greatly 
increases ~,. Filaments containing different amounts of thoria but 
otherwise alike seem in general to give values of », proportional to the 
amount of thoria. 

Amounts of thorium metal present in the filaments. We see from Eq. 
(29) in the case of parabolic distribution, that the average concentration 
throughout the filament is one-half that at the center and thus according 
to Eq. (21) it is equal to GR/4. In the stationary state we may replace 
G by G, which we find from Eq. (35), and thus taking R=0.00389 cm 
and by converting from atoms per cm® to parts by weight by dividing 
by 5.010” we obtain the following average concentrations of metallic 
thorium in parts per million. 

1800° 1900° 2000° 2100° 2200° 2400° 2800° 3000°K 

1.4 2.8 5.0 8.4 14 32 118 200 

Relation between induced evaporation and 6. The activation curves 
obtained in Exp: 528 after the filament had been well aged gave straight 
lines when 6,, —@ was plotted on semi-logarithmic paper against the time, 
i.e., the data agreed with Eq. (10) which was derived on the assumption 
that G is constant. The foregoing analysis, however, indicates that in 
the neighborhood of 2100°, G should change considerably during the 
run if G is initially very different from the stationary value G,. 

For example in an activation run at 2100° (first part of run U,) which 
followed after a flashing at 2800° for 30 seconds, 6 increased from 0.01 
to 0.83 in 20 minutes. The line obtained by plotting @..—6 on semi- 
logarithmic paper was so nearly straight that its slope did not change 
more than 2 per cent in the course of its length, while calculation by 
Eq. (23) shows that G should have decreased about 6 per cent. Run 
N, furnishes a more striking example for the slope remained constant 
during the whole run of 277 minutes at 2055° although G decreased 55 
per cent during this time. 

This result may be explained by assuming either that @ is not strictly 
a linear function of log i, or that the induced evaporation is not quite 
proportional to 6. Taking the latter assumption as the more probable, 
it has been found that the data of runs U, and M,, as well as several 
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others of a similar nature, are very satisfactorily accounted for if the 
factor (1—@) in Eq. (8) is replaced by 

f(0) =1—0.82 0—0.18 6°, (38) 
the last 2 terms representing the induced evaporation. 

Thus when @ is small, the chance is 0.82 6 that an atom of thorium 
in arriving at the surface shall cause the evaporation of one already on 
the surface. But when @ is nearly unity, the chance that an atom of 
thorium arriving at the surface can find a place in the surface without 
driving off one already present is 1.36 (1—@). This suggests that the 
thorium atoms are not closely packed in the surface film even when 
é=1. 

Normal evaporation of thorium. ‘The actual rate of evaporation of 
thorium may now be calculated from the experimental data by Eq. (8). 
For the present, however, we do not need to assume that this rate is a 
linear function of @ and will therefore write the equation in the form 

No dé/dt=DG f(6)—-E, (39) 
where f(@) is given by Eq. (38) but with sufficient accuracy for most 
purposes may be placed equal to (1—6). 

As an illustration of the method used in calculating the values of E 

as a function of 6, let us consider the data of run X, given in Table V. 


TABLE V 
Rate of evaporation of thorium as a function of @ 


Data from Run X,, Exp. 528 Filament temperature 2300°K 
D=1.30X10-° Gr = 1.09 X 10! R=0.00389 cm 


(1) (2) 
tm 6 








(3) (4) (5) (6) 
104d /dt f(6) 10-“DG 10-“E 
(atoms cm-* (atoms cm-? 
(min.), (per sec.) sec."!) 


0 7.86 
2 9.07 
+ 9.60 
9.98 


6 

10.36 
10.58 
10.81 
11.26 





8 

10 .< 
12 . 296 
16 .201 
26 .141 
24 -111 
28 .094 


11.64 
11.94 
12.25 


cos s9s5 9999 








The values of @ were calculated by Eq. (1) from the saturation currents 
i shown in Fig. 4. The derivative d@/dt was obtained much more accu- 
rately than would be possible from a direct plot such as Fig. 4, by plotting 
(6—@..) on semi-logarithmic paper against ¢ so as to obtain a curve which 
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was very nearly straight, whose slope could thus be accurately measured 
at any point. According to Eq. (10) this slope is equal to k which is 
here considered variable and it is evident that 

d0/dt = —k(0—6..) (40) 
The value of the derivative found in this way is independent of the 
choice of @,,. 

The values of DG in Col. 5 are calculated by means of Curve III in 
Fig. 6, in the following manner. In run V, the filament had been heated 
for several hours at 2250° until G had reached its stationary value 
G, which by Eq. (35) is equal to 8.810” at 2250°. In run W, which 
followed immediately after V;, the filament was at 2105° for 196 minutes 
during which G fell from 8.8 x 107° to 6.2 x 10° while @ rose from 0.12 to 
0.857. The filament was then maintained at 2000° for 47 minutes which 
lowered G to 6.05 X 107° and raised @ to 0.883, and this was the condition 
of the filament at the beginning of run X, which followed. The initial 
value of DG at the head of Col. 5 (Table V) was obtained by multiplying 
G by 1.30X10-*, which is the value of D at 2300° from Table IV. 

From ¢t», the time in minutes in Col. 1, we obtain Dt/R® and from 
Curve III in Fig. 6 find the corresponding values of RG/C,. These are 
proportional to DG—DG,, and thus, since we know the initial value of 
DG and obtain G,=1.09X10* at 2300° from Eq. (35), we calculate the 
remaining values of DG in Col. 5. 

The rates of evaporation E given in the last column are calculated by 
Eq. (39) from the corresponding values of d@/dt, DG and f(@) in columns 
3, 5 and 4; No being given by Eq. (15). Although @ decreases in the 
ratio 10 to 1, E decreases by only 15 per cent, and over most of the 
range may be taken as constant. 

In run S;, which is also shown in Fig. 4, a similar analysis indicates 
that DG at t=5 minutes, was 5210" and fell to 15X10" at ¢=106. 
Notwithstanding these much higher values of DG, the values of E were 
not essentially different, averaging 16.310" in the range 6=0.7 to 
0.2 while a similar average in run X; gives 11.5 X10". This difference 
is equivalent to a change in temperature of 18° at 2300°. In run S; 
also, E was practically independent of 0. 

The evaporation data from run X, are probably much more reliable 
than those from run S;, since in the former run the term con- 
taining DG in Eq. (39) was much smaller in comparison with the 
do/dt term than in the latter run. The ratios of the values of these 
terms at the beginnings of these runs were 1:10 and 4.5:1 respectively. 
Of course the factor DG during the course of a run is subject to con- 
siderable uncertainty as it is not determined directly but is calculated 
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from the data of Fig. 6 and involves assumptions as to the distribution 
of thorium throughout the filament. The term involving d@/dt on the 
other hand is found from the experimental data. Therefore it is felt 
that very little weight should be given to determinations of E in runs 
such as S; in which G had high values. 

Two other methods were used in some experiments to determine values 
of E. These may be illustrated by the data of run V;. The filament 
was first heated at 2800° for 30 seconds, and then run at 2195° for 60 
minutes, at which time the temperature was changed to 2250°. Between 
the first and the 12th minute at 2195°, @ increased from 0.23 to 0.88, 
and the values of (@—@,.) gave on semi-logarithmic paper a straight 
line whose slope corresponded to k =35.5X10~ sec.-!. The value of 6,, 
which had to be chosen to get this straight line was 0.943. Substituting 
these values in Eq. (14) gives »+v=1.53X10". Assuming, in accord 
with the result of run X, that »=0 in comparison with », Eq. (13) 
becomes DG = Nok, which gives DG = 26.8 X10", or taking D =4.85 x 10-!° 
at 2195° by Eq. (34) we find G=5.53 X10", and E=1.53 X10". 

During the remainder of the run at 2195°, 6 increased to a maximum of 
0.921 at ¢=24 minutes and then decreased slowly as in run S;, Fig. 4. 
At t=60 min. the data showed 6=0.882 and d6é/dt = —0.20X10- sec.-!, 
the filament temperature still being 2195°. When at this point the tem- 
perature was changed to 2250° the filament was rapidly deactivated, 
so that the curve of @ as a function of ¢ gave a sharp kink, the slope 
changing abruptly from —0.2010~* to —4.5X10~*. Substituting these 
data obtained just before (1) and just after (2) the change of temperature 
into Eq. (39) together with the value of No from Eq. (15) we get the 
equations, 

—0.15 X10" = D,Gf(6)—E, 
and —3.40 X10" = D.Gf(6) — E, 
where D,=4.85X10-', and D,.=8.22X10", and f(@)=0.153 (for 6= 
0.882). Solving these equations we find 
G =1.35 X10'°E; —0.20 x 107 
and E,=3.15 X10"+1.69E; - 

In this way by measuring the change in slope caused by a change in 
temperature, the value of E at one temperature can be found from that 
at another. If in this example we take the value E,;=1.53 X10" found 
from the first part of run V; we obtain E,=5.74X10" as the rate of 
evaporation at 2250° and G=1.86 X10 which is about one third the 
concentration gradient found at the beginning of the run. 

Table VI contains a summary of the most reliable data on the rate of 
evaporation in Exp. 528. Col. 2 indicates the method used to determine 
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E, the three methods which have been described above being designated 
respectively I, II and III. The letter A or D is used to indicate whether 
the filament became activated or deactivated during the run. Col.%4 
describes the condition into which the filament has been brought by its 
previous treatment; it gives the value of G observed in the middle of the 
range of @ covered by the run. In some of the runs @ varied over a wider 
range than indicated in Col. 5, but the value of E given in the last 
column is based on the observations in the range given and in this range 
E seemed to be substantially independent of @. 


TABLE VI 


Summary of data on rate of evaporation of thorium. Experiment 528 





Method Temp. 10-2G Range of 6 


I—A 2055°K 
II—A 2105 
I-—D 2150 
I-—D 2195 
II—A 2195 
IlI—D 2250 
I— 2300 
I— 2300 
2400 
2450 
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In Fig. 7 the logarithms of E are plotted against the reciprocal of the 
temperature. It is seen that the points lie along a straight line which 
corresponds to! 


logio E = 31.434 —44500/T. (41) 
Values of E calculated from this equation at intervals of 100° are given 
in Table IV. The slope of the line in Fig. 7 is 44500 degrees which cor- 
responds to a heat of evaporation of 204000 calories per g-atom. 
Although in the range of 6 covered by the foregoing runs, E is found to 
be nearly constant at any given temperature it is certain from kinetic 
considerations that E must become proportional to @ for sufficiently 
small values of @ where each thorium atom must act independently of 
others. In fact a few runs which had been carried to particularly low 
values of 6 gave some indication of a falling offin E. In run U, at 2450° 
E was practically constant between 0.2 and 0.8 but at @=0.1 had fallen 
to 65 per cent and at @=0.07 to 48 per cent of its value at @=0.2. 


“In the abstract in the Phys. Rev. 20, 107 (1922), the first term in the second 
member of this equation was given as 7.76. In the early work » was used to express 
evaporation in terms of gram atoms Xcm- sec.“ In the abstract the definition of » 
was given in atoms instead of gram atoms but by an oversight the equation was not 
correspondingly altered by multiplying by the Avogadro constant 6.1 x 10”. 
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To get more information as to the evaporation in this region of low 6 
about 30 runs were carried out in Exp. 528 (after run X,) with values of 
6 ranging from 0.02 to 0.3. The data were very consistent but in view 
of more recent experience with other types of active cathodes, it seems 
probable that they were subject to considerable errors due to incomplete 
evaporation of thorium from portions of the filament near the leads which 
were at a somewhat lower temperature. It is evident, because of the 
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Fig. 7. Evaporation and diffusion of thorium as functions of temperature. 


enormous ratio between the emission from active and deactivated fila- 
ments, that if a very short section of the filament retains its thorium, 
serious errors may arise in the determination of @ particularly for low 
values of @. With filaments in a fairly active condition low testing 
temperatures such as 1450° can conveniently be used. Since the cooling 
effect of the leads is much more pronounced at these low temperatures 
than at the higher temperatures at which deactivation occurs, there 
is no appreciable emission from end portions of the filament which may 
have retained their thorium. At the higher testing temperatures used 
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with the low values of 6, more serious errors may occur. In future ex- 
periments all these difficulties can be avoided by dividing the anode into 
sections on the guard-ring principle.” 

The general results of these runs with low values of @ are as follows. 
At temperatures from 2200° to 2250° E did not seem to fall off in value 
as @ decreased from 0.3 to 0.07 but showed a slight tendency to rise. 
At higher temperatures E did decrease with 0, the more rapidly the higher 
the temperature, but even at 2400° decreased less rapidly than in pro- 
portion to @. As a result of these variations of E the temperature 
coefficient of EZ, and the heat of evaporation, are less than at higher 
values of @. 

The activation runs such as those of Tables I and II furnish some 
data on evaporation at high values of 6. Assuming E to be a linear 
function of 6, we see by Eq. (7) that the value of E for @=1 is m+», 
and this may be denoted by Z£;. Then Eq. (14) takes the form 

1—6,, = E;/Nok (42) 

Therefore (1—8@,,) plotted against 1/k should give a straight line pass- 
ing through the origin. The observed values of 6, in Table II fulfil 
this condition satisfactorily. The slope of the line gives E,/No from which 
we find E,=1.66X10". The 6th Column of Table II contains values 
of 6. calculated by Eq. (42) with this value of EZ, and the data for k 
in Col. 4. The agreement with the observed values is probably within 
the experimental error and justifies the assumption of a linear relation 
between E and @ in the region from @=0.94 to 1.0. The value of E, 
from the data of Table II is 2.75 times greater than the value of E 
calculated from Eq. (41) for this temperature (2055°). 

The data in Columns 4 and 5 of Table I were also used to calculate 
values of E, by Eq. (42). These were compared with E calculated for 
the same temperatures from Eq. (41) and were found on the average to 
be 3.12 times greater. Col. 6 Table I contains values of 0, calculated 
by Eq. (42) from the data for k in Col. 4 by placing E,=3.12E where E 
was calculated by Eq. (41) for the temperature 7, (Col. 2). The agree- 
ment between the observed and calculated values of 6,, is excellent. This 
is of especial significance when it is remembered that the data of Table I 
were obtained with a filament containing 1.5 per cent of thoria while the 
filament of Table II had 1.0 per cent. 

The work at low, intermediate, and high values of @ thus leads to the 
following conclusions. In the range from 6=0 to about 0.1, E increases 
probably about in proportion to 6, but before 0.2 is reached increases 
much more slowly. Between 0.2 and 0.8, E is approximately constant, 


% Langmuir and Kingdon, Science, 57, 58 (1923) 
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showing however a tendency to increase slightly with 6 at temperatures 
as high as 2400°, but to decrease slightly as 6 increases at temperatures as 
low as 2200°. Above 0.8, E increases again and between 0.85 and 1.0 
it increases linearly with @ until at @=1, E is about 3 times as great as 
in the range from 0.2 to 0.8. 

Variation of heat of evaporation with 6. Dushman"™ has proposed the 
following empirical equation for the velocity of monomolecular reactions 


éhse Cpe WAT (43) 
dt 

where c is the concentration of the reacting substance, and hy is the 
heat of reaction per molecule. Although this equation is certainly not 
rigorous it has been found to apply with a satisfactory degree of approxi- 
mation to many fundamental processes which depend on temperature. 
The first member of the equation is the reciprocal of the average life r 
of the molecules of reacting substance, and in this form the equation 
becomes applicable to many classes of phenomena for which we have 

previously had no similar relation. 
The average life of a thorium atom on a surface of tungsten is given by 
t=N0/E (44) 
Combining these equations and substituting Q= Nh»v, where Q is the 
heat of evaporation of thorium from tungsten per gram atom, we obtain 
E=(NoQ0/Nh)e®/®? (45) 
When the values of No, N and h are substituted in this equation it 

becomes 

E=1.804 X10'7Q6e-°/* (atoms sec.~! cm~*). (46) 
We may now calculate the value of Q as a function of @ from our data 
on the rate of evaporation as illustrated in Table VII. The second 
column gives the rate of evaporation at 2300° calculated from Eq. (41) 
which merely summarizes the experimental data. The value at @=1 
is made 3 times that in the range 0.2 to 0.8 in accordance with our 


TABLE VII 
Heat of evaporation of thorium as a function of 0; filament temperature 2300°K 








Rate of evaporation Heat of evaporation 
E 






(atoms sec.-! cm-*) (calories per g-atom) 








0.2 1.2210" 182,800 
0.5 1.22 « 187,100 
0.8 1.22 « . 189,200 
1.0 3.66 “ 185,200 








%S, Dushman, J. Amer. Chem. Soc., 43, 397 (1921) 
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previous conclusion. The 3rd Col. gives values of Q calculated from 
these data by Eq. (46) but converted to calories per gram atom. 

Considering the approximate nature of the Dushman equation the 
agreement between the calculated values of Q and the value Q = 204,000 
found more directly from Eq. (41) is satisfactory. It is probable that the 
differences in Q for different values of @ are accurate to much higher de- 
gree. Almost the same differences may be calculated from the same data 
by the Boltzmann equation. The most noteworthy feature of these 
variations of Q with @ is that they are so small. It indicates that the 
forces acting between adjacent thorium atoms on the surface is very 
small compared to those involved in holding these to the underlying 
tungsten atoms. The increase in Q between 0.2 and 0.8 corresponds to 
attractive forces while the decrease between 0.8 and 1.0 is evidence that 
the atoms are subjected to repulsive forces by the neighboring thorium 
atoms when they become more crowded. The effects are so small, how- 
ever, that they prove that the thorium atoms are not closely packed. 
The results are thus in accord with the conclusion we have reached in 
regard to the magnitude of Np. 

Volmer and Estermann" have shown that mercury atoms that con- 
dense on a solid mercury crystal may retain their mobility of thermal agi- 
tation in an adsorbed film, as a kind of 2-dimensional gas. It is highly 
probable that thorium atoms on a tungsten surface at sufficiently high 
* temperature are able in a similar way to move over the surface. The 
phenomena of induced evaporation, however, seem to indicate that this 
degree of mobility is very limited, for a thorium atom coming from 
within the filament to a surface only partly covered with thorium, would 
otherwise be able to push a thorium atom on the surface into an adjacent 
position instead of inducing it toevaporate. The average life of a thorium 
atom in a given position on the tungsten surface lattice (with respect to 
surface mobility) must therefore be long compared to the life of a thorium 
atom (with respect to evaporation) which has a thorium atom under it. 

A summary of the properties of thoriated tungsten filaments is con- 
tained in Table VIII. With any such filament we may classify the 
phenomena into 2 groups: (I) those characteristic of the stationary 
state in which the metallic thorium is diffusing out of and leaving the 
filament at the same rate as it is produced within the filament by reduc- 
tion of thoria, and (II) transient phenomena due to departures from the 
stationary state caused by changes in temperature or removal of thorium 
from the surface of the filament by sudden flashing or positive ion 


17 Volmer and Estermann, Zeitschr. f. Physik, 7, 1 and 13, (1921) 
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bombardment. These transient phenomena are of two kinds: (Ila) 
Those in which the rate of production of thorium metal in the filament 
is not equal to the rate of diffusion to the surface (so called steady state), 
and (IIb) Those in which the rate of evaporation from the surface is 
not.equal to the rate of arrival at the surface. 


















TABLE VIII 
Properties of a thoriated tungsten filament in the steady state 
Filament contains 1 per cent of ThO, Diameter 0.00389 cm 














(1) (2) (3) (4) (S) (6) (7) 
Electron Time of Time of 
i Gr 6r emission Life recovery activation 
tr T 1/a 1 
(atoms cm-*‘) (amp. cm-*) (hours) (sec.) (sec.) 

1300°K 6.61X10'7 0.99997 4,14x10-4 1.50 101° 
1400 2.2310 0.99975 3.12x10-3 1.11109 1.4310"! 
1500 6.4110" 0.99878 0.0179 1.17 x 108 5.23 X 10° 
1600 1.6110 0.99528 0.0812 1.63 < 107 2.88 X 108 
1700 3.6410! 0.9848 0.287 2.86 X 108 2.15107 
1800 7.5110! 0.9605 0.772 720,000 6.10 10° 2.02 x 108 
1900 1.43107 0.9191 1.59 94,000 1.53105 2.18 X 10° 
2000 2.57107 0.8713 2.89 15,100 44,100 2.48 X 104 
2100 4.36102 = 0.781 3.43 2,897 14,300 9.48 x 10 
2200 7.03107 0.551 1.24 643 5,130 2.10 108 
2300 1.091022, 0.139 0.114 164 2,010 531. 
2400 1.63102 0.0601 0.168 47 855 7.0 
2500 2.3610, 0.0355 0.357 14.6 389 1.28 
2600 3.31102 0.0207 0.774 5.01 191 0.268 
2800 6.0810 0.0088 3.48 0.74 51 0.0161 
3000 1.031072 0.0041 13.5 0.14 17 0.0014 























Col. 2 in the table gives values of G, (the concentration gradient at 
the surface when the stationary state (I) has been reached) calculated 
by Eq. (35). Col. 3 gives 6,, the fraction of the surface covered with 
thorium when the filament is in the stationary state (I). These data 
are calculated by means of Eq. (12), 6. being replaced by 6, and G by G,. 
The values of D are taken from Table IV or Eq. (34). In accord with 
the experimental data discussed previously, the following convention 
has been adopted for values of vp and ». 


Range of @ Yo v 

0.0 to 0.1 0 10 E 
0.1 to 0.85 E 0 
0.85 to 1.0 —10.33 E +13.33 E 


The value of EZ is given by Table IV or Eq. (41). 
Col. 4 gives the electron emission 7, in amperes per cm? at the tempera- 
ture 7 (Col. 1) when the fraction 6, of the surface is covered with thorium, 
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Dushman* has shown that the electron emission from any heated metal 
is given by 

i=60.27%e*/? amp. cm~, (47) 
where } is a constant which depends on the material of the cathode 
or the condition of its surface. 

Eq. (2) may be written 

b=bo+0(b1 — bo) (48) 
Dushman has obtained the values: bo = 52,600, b: = 34,100. Substituting 
these in Eq. (48) and then placing the value of 6 in Eq. (47) we obtain 
the data given in Col. 4. 

The “‘stationary’”’ emission from a thoriated filament thus reaches a 
maximum somewhat above 2000° and then falls gradually to that of 
pure tungsten. 

Col. 5 gives the life of the thoria in the filament calculated from 
pb, (Eq. (37)) by the relation r=1.22 X10""/p,.. This represents the time 
required for the total thoria content of the filament to fall to one eth 
of its original value. 


Col. 6, under the heading Time of recovery, gives the time required 
for transient effects of the first class (IIa) to fall to one eth of their value. 
Thus if the filament has been kept at 2000° until the stationary state has 
been reached, and the temperature is then changed to 2100°, G begins 


to change from the value 2.57 X 10° and approaches the new stationary 
value 4.36 X 10°, in such a manner that it takes 2897 seconds to go to 
within one eth of the way to the final value. The quantity a, whose 
reciprocal is the time of recovery, is calculated and defined by Eqs. 
(23) and (24). The time of recovery gives a measure of the time 
needed to bring the filament into the stationary state. 

The last column of the table gives the time of activation, the reciprocal 
of the rate of activation k, defined by Eq. (10). The values in the table 
were calculated by means of Eq. (11), replacing G by G,. This time of 
activation is the time needed for @ or 7 to return to within one eth of the 
way back to the “steady” value (@,, or i.) after any small change pro- 
duced by altering the temperature for such a short time that G, is not 
appreciably affected. The time of activation is also a measure of the 
sensitiveness of the filament to effects which tend to destroy the thorium 
film on the surface, such as positive ion bombardment or chemical 
action of residual gases. 

The data of Table VIII apply strictly only to a filament of the diameter 
given. By Eq. (36) we see that G, should increase in proportion to R 


18 Dushman, Phys. Rev.£20, 109 (1922), and Gen. Electric Rev. 26, 154 (1923) 
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and this will change 6, and i,. The life 7 is independent of R. By Eq. 
(24) it appears that the time of recovery is proportional to R’, and by 
Eq. (11) the time of activation is inversely proportional to R at low 
temperatures and independent of R at high temperatures. 

By introducing carbon into the filament, as for example by heating 
the filament in hydrocarbon vapors at low pressure, the value of G, 
may be greatly increased, causing a corresponding increase in 6, and i, 
but a decrease in the life and in the time of activation. The decreased 
sensitiveness to gas effects may be of very great advantage in practical 
devices utilizing this type of cathode. 

The writer wishes at this time to express his appreciation of the as- 
sistance and suggestions he has received from Dr. William Rogers, Mr. 
S. Sweetser and Dr. S. Dushman. 
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DIELECTRIC LOSSES AT RADIO FREQUENCIES IN 
LIQUID DIELECTRICS 


By A. B. Bryan 







ABSTRACT 







Dielectric losses and dielectric constants at radio frequencies for nitro- 
benzene, water and xylene.—The method of resistance variation was used to 
measure the phase difference y and dielectric constant K for frequencies 
between 210° and 14X10° cycles/sec. Special cells were required. (1) 
Variation with frequency. The results agree approximately with the equations: 
For carefully dried nitrobenzene at 30°C, y =.028°+6.03 x 10*/f; for distilled 
water at 23.5°, y =0.8°+2.09 x 10°/f. These indicate that in addition to the 
true dielectric loss there is a leakage through the liquid proportional to 1/f. 
For xylene, ¥ was too small to measure, less than .01° at 310° cycles. K 
was found to be practically independent of the frequency, being 2.24 for 
xylene and of the order of 100 for water. (2) Variation with temperature, for 
nitrobenzene. K decreased from 42 at 20° to 24 at 14.2°, while y increased in 
the same range in the ratio of 7 to 1. These values were obtained, however, 
for a sample of nitrobenzene for which y was 12 times as great as for a carefully 
dried sample. 
















INTRODUCTION 









HEN an alternating electromotive force is applied to a condenser 
a certain amount of power loss occurs due to leakage through 
or over the surface of the dielectric and to hysteresis, or inter-molecular 
friction, or some such phenomenon, in the dielectric material. This 
latter source of power loss is usually denoted by the term dielectric loss. 
Power loss in a condenser is conveniently expressed in terms of phase 
difference which for small values may be obtained from the relation, 
¥=RwC 
where R and C are the equivalent series resistance and capacity, re- 
spectively, of the condenser and w is 2x times the frequency. For 
small power losses the phase difference of a condenser is the same as its 
power factor. 

A large number of measurements of phase difference have been made 
on solid dielectrics of various kinds. Reference may be made to the 
work of Schott! and to the more recent investigations of the Bureau 
of Standards’ and of MacLeod.’ In the present work the measurements 
are extended to liquid dielectrics. 




















1 Schott, Zeit. fur drahtlose Telegraphie und Telephonie. 18, p. 82 
* Technological Papers of the Bureau of Standards, No. 216 
3 MacLeod, Phys. Rev. 21, 53, Jan. 1923 











































A. B. BRYAN 


EXPERIMENTAL METHOD 


The experimental method is one which has usually been used for 
previous measurements on solid dielectrics. A diagram of the circuit 
is shown in Fig. 1. The measuring circuit LCT is inductively coupled 
to L’, which is part of a radio frequency oscillator. Heavy movable 
copper links were provided by which either the variable air condenser C, 
or the condenser to be tested C could be connected in the circuit. The 
high frequency current in the circuit was measured by means of the 
thermocouple T and the micro-ammeter G. Mercury cups MM were 
provided which could be connected by special high frequency resistance 
units of any desired value, thus making it possible to vary the resistance 
of the circuit. Each of these resistance units was made of a short straight 
piece of manganin wire, suitably mounted and of sufficiently small 
diameter to make negligible the change in resistance with frequency. 
The complete circuit was supported on ebonite posts several inches from 
the table to decrease losses due to leakage and to capacity to ground of 








Fig. 1 





various parts of the circuit, and was mounted on a rolling platform 
so that its distance from the oscillator could easily be changed. By 
thus varying the coupling between the oscillator and measuring circuit 
any convenient current could be obtained in the measuring circuit. 
The condenser C, is a General Radio Company ‘Precision Condenser’’ 
with micrometer screw adjustment and a maximum capacity of 1500 
micro-micro-farads. Condensers of various forms were used for the 
condenser C; they will be described when the results are considered. 
All measurements of frequency were made with a General Radio Com- 
pany ‘Precision Wavemeter.” 

In making a measurement, C is first connected in and the oscillator 
tuned to resonance with the measuring circuit. The total resistance of 
the measuring circuit, including the equivalent series resistance of C, 
is then determined by resistance variation. C is next replaced by C; 
which is varied until oscillator and measuring circuit are again in re- 
sonance and then the total circuit resistance is again determined as 
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before. The difference between this second value of the measuring 
circuit resistance and the value first obtained is the equivalent series 
resistance R of the test condenser C, the equivalent series resistance of C; 
being negligible. The capacity of C is obtained from the scale reading 
of the calibrated condenser C;. The frequency having been measured 
with the wavemeter, the phase difference may then be calculated. 


RESULTS 


Some preliminary measurements were made on several condensers 
having solid dielectrics. Values of phase difference and dielectric 
constant were obtained for frequencies between 200,000 and 1,200,000 
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cycles per second and both phase difference and dielectric constant were 
found to be approximately constant over that frequency range, in 
general agreement with the results of previous investigators. 

The liquid dielectrics investigated were nitrobenzene and distilled 
water. Figs. 2 and 3 show the variation of y and the dielectric constant 
K with temperature and with frequency for nitrobenzene. y is expressed 
in degrees. The condenser used in these cases is shown in Fig. 4. The 
two plates were concentric silver cylinders A and B fixed in a glass cell. 
The cell had a ground glass stopper to prevent the nitrobenzene from 
absorbing moisture from the air. The cylinders had diameters of 3 
and 4 cm and were 5 cm long. Variations of temperature were obtained 
by dipping the cell in an oil bath which was electrically heated. The 
absolute values of K for nitrobenzene are only approximate since the 
dimensions of the cylinders were not accurately known. 
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Considerable trouble was experienced in removing all the moisture 
from the nitrobenzene. The method first tried was to let the nitrobenzene 
stand some time over calcium chloride and then filter and distill it. 
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The method finally used was to add red hot anhydrous copper sulphate 
and then distill. The nitrobenzene used to obtain the results shown in 
Fig. 2 had not been as thoroughly dried as that which was used for the 
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results in Fig. 3, hence ¥ at 30°C and 625,000 cycles is 1.6° in Fig. 2 


and only 0.13° in Fig. 3. It is possible that if the nitrobenzene were 
perfectly dry the values of y would be still smaller. 
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The results obtained for distilled water are given in Fig. 6. The 
type of condenser which was finally adopted for distilled water is shown 
in Fig. 5. Two square platinum plates A and B, 1.6 cm square and 0.3 
cm apart, formed the plates of the condenser. They were supported by 
platinum wires which ran down through the ground glass stopper of the 
glass cell. The plates dipped into the distilled water which was 
in a small platinum crucible C in the bottom of the cell. Several other 
condensers were first tried out, including the one which had been used 
for nitrobenzene, but it was found that the conductivity of the water was 
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Fig. 6 
too large until the above condenser was used, where the water comes in 
contact with platinum only. The distilled water used was obtained from 
an electrically heated still with a tin condensing coil. 

In Fig. 6 the values of the total capacity C instead of the dielectric 
constant K are plotted. It is not possible to calculate K accurately 
because the dimensions of the plates cannot be measured accurately 
and because of the difficulty of allowing for the capacity of the backs of 
the plates and of the connecting wires. But if we assume that the given 
dimensions are correct and that the effect of the backs of the plates and 
of the connecting wires is to increase the capacity by 50 per cent, then 
we get K=97 approximately. At any rate, water has a very large di- 
electric constant at the frequencies used. 

Several measurements were also made on xylene with another type 
of condenser but the losses in this liquid were found to be too small 
to be measured with the method used. For xylene the phase difference 
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is less than 0.01 degree at 300,000 cycles. K was independent of the 
frequency and equal to about 2.24. 

The fact that y decreases rapidly with increasing frequency in the 
case of both liquid dielectrics, indicates that a large part of the power’ 
loss is due to leakage through the liquid. If this leakage is such as 
would be produced by a constant parallel resistance R, then it may be 
shown that 

y=1/22R,fC 
and the y due to this source will be inversely proportional to the fre- 
quency f. If we assume that y also has a constant component yo due 
to true dielectric loss, then the total phase difference will be 

Y=YWtA /f. 
The curves for y in Figs. 3 and 6 are plotted from this equation, the 
constants being of course determined by selecting two points through 
which the curve shall pass. For water y=0.8° and A=2.09X10°. 
For nitrobenzene ¥=0.0282° and A=6.03X10*. The agreement of 
the experimental points with the curve is good in the case of water and 
only fair for nitrobenzene. In the latter case the relation between y 
and f might easily be of the form 

y=A-—Bf. 
It therefore seems probable that in liquid dielectrics, as well as in solid 
dielectrics, the phase difference due to true dielectric loss and also the 
dielectric constant are practically independent of the frequency. 

In conclusion, the writer wishes to express his gratitude to Professor 
H. A. Wilson, under whose direction this work was done. 
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POLARIZATION EFFECTS SHOWN BY FILMS OF 
CERTAIN FUSED SALTS 


By D. B. DEopHAR AND G. B. DEODHAR 


ABSTRACT 


Some experiments made by Talbot and Lord Rayleigh have been repeated 
and extended. A film was prepared by evaporating a solution to dryness on 
glass, then fusing the salt and letting it cool. In the case of potassium nitrate 
and potassium chlorate, when the film was placed under a microscope and 
between crossed nicols, Talbot and Rayleigh had observed irregularly dis- 
tributed patches varying from dark to light, and variously colored. It is now 
found that silver nitrate and benzoic acid also show the effect. The intensity 
diminishes with time but the variations can still be observed after six months. 
A film formed by evaporating a binary solution of a mixture of any pair of the 
following: alum, boric acid, borax and citric acid, also showed the effect al- 
though films containing only one constituent were inactive. The effect is 
doubtless the result of mechanical strains developed on cooling the film, similar 
to those obtained when glass is chilled from above its softening point. 


INTRODUCTION 


OME years ago the late Lord Rayleigh! had published a note upon a 
remarkable phenomenon of Newtonian colors seen by Talbot, when 
fused nitre was observed under crossed nicols by means of a polarizing 
microscope. Talbot found that the fused film appeared to be made up 
of crystalline plates of irregular shape, often fitted together like a tes- 
selated pavement. If one plate was darkened by rotation of the nicol, 
others remained visible in varying degrees of brightness. On suitably 
increasing the quantity of the salt he could see the display of gorgeous 
colors. He also observed that the whole appearance rapidly changed in 
a surprising manner, and that the whole color phenomenon ceased at the 
end of about a quarter of an hour, the field becoming quite dark. The 
late Lord Rayleigh repeated this experiment with nitre and extended his 
observations upon potassium chlorate in practically the same way as 
Talbot had done, and found that the chlorate colors were less beautiful 
than those of the nitre. He also found that the potassium chlorate colors 
lasted for a longer time than nitre colors. Reinders,? Brauns,* Lehmann‘ 


1 Rayleigh, Collected papers, Vol. 6 

* Reinders, Zeits. phys. Chem. 1912 

* Groth, Physikalische Krystallographie 
Lehmann, Molekular Physik 
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and Tutton® have also touched upon these color phenomena. In the 
present work the observations of Talbot and Rayleigh have been repeated 
and extended. 


OBSERVATIONS 


(a) A few drops of potassium nitrate solution were laid upon a clean 
piece of glass, and were slowly evaporated to dryness. Then the dry 
substance was melted, the glass piece was removed and allowed to cool 
for a little while, and after some four or five minutes it was placed 
on the stage of a microscope having a nicol polarizer below and a nicol 
analyzer within the body of the instrument, set for extinction. Observa- 
tions were made with sunlight and the whole field presented the appear- 
ance of a tesselated pavement of irregular boundary lines as described by 
Talbot. Different patches, some dark, some bright, some of varying 
shades of brightness and darkness, and some exhibiting brilliant New- 
tonian tints of varying intensity, were irregularly distributed through- 
out the field. On rotating the polarizing nicol through 90°, it was 
observed that the dark patches became bright, while the color tints 
changed to complementary ones. As time increased, the color effect 
appeared to decrease; but its rate of decrease was found to be slower and 
slower as time went on. After the lapse of about an hour’s duration, 
changes in the intensity of the colors could not be detected. 

(b) Potassium chlorate was treated in the same manner and its general 
behaviour was found to be the same as that of potash nitrate. The 
only difference detected was that the intensity of color was far stronger 
than that with nitre, contrary to the result obtained by Rayleigh. The 
rate of decrease in the color intensity with increase in time was found to be 
less than in the case of nitre. 

(c) Silver nitrate was tried and found to behave as above. 

(d) Benzoic acid showed very gorgeous colors under crossed nicols, 
which changed very slowly with time. It was further noted that if the 
heat given is just sufficient to melt the acid, the colors are very vivid; 
but if the heating is pushed a little too far, the color effect becomes quite 
dull. 

(e) Another peculiarity which is observed during these experiments is, 
that if while observing through the eyepiece, under crossed nicols, a light 
scratch is made upon the fused substance on the glass piece by means of a 
pin point, at once colors begin to jump out of the point of disturbance in 
all directions, the small portion round the scratch giving an appearance 
of a peacock feather. The path of the pin point could be followed by the 


5 Tutton, Text book on crystallography 
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beautiful radiant color streaks. This effect is most striking a few minutes 
after melting; and when the substance becomes cold, the pin point cannot 
scratch and consequently the radiant color streaks cannot be observed. 

(f) In order to find out the effect of time on our substances, they 
were carefully preserved. Unlike Talbot and Rayleigh, we found that 
all our preparations showed colors even at the end of six months. Of 
course, the intensity was not the same as that observed when the speci- 
mens were freshly prepared, but the effect was still clearly evident. 
This shows that the rate of diminution of intensity becomes extremely 
slow as time advances. These specimens are still being watched. 

(g) A few drops of a solution of ordinary alum in water on clean glass 
were slowly evaporated to dryness, and the behaviour of the substance 
was observed without fusing it; but nothing could be seen under crossed 
nicols. Citric acid, boric acid and borax also showed the same negative 
behavior. However in the case of a mixture of boric acid and alum, 
color effects were observed, which showed the usual complementary 
relations on rotating the polarizing nicol. A mixture of citric acid and 
boric acid and of borax and citric acid also showed the same comple- 
mentary color relations. 


DISCUSSION 


Our observations of color effects indicate the existence of anisotropic 
molecular aggregates with their optic axes oriented at random; and it 
appears that the double refraction property may be due to the inter- 
molecular strains and stresses developed on account of the setting of two 
different types of molecules in mixtures; similarly as a result of the cool- 


ing of the fused substances, unequal stresses and strains would be in- 
troduced. 
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